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Abstract 

The  tracking  performance  of  the  Fire  Fly  manual  director 
gunsight  is  evaluated  during  the  terminal  phase  of  air-to- 
air  gunnery  using  the  F-106  aircraft  as  the  test  bed.  A 
digital  simulation  of  the  closed  loop  tracking  task  is 
performed. using  the  EASY  Dynamic  Analysis  programs  developed 
by  Boeing  Computer  Services,  Inc. ^ The  multi-axis  analytical 
pilot  model^developed  for  the  Terminal  Aerial  Weapon  Delivery 
Simulation  (TAWDS)  program  by  McDonnell  Aircraft  Company,"^ is 
adapted  for  use  in  the  F-106/Fire  Fly  gunsight  simulation. 

A set  of  nominal  pilot  model  parameter  values,  determined  by 
root  locus  analysis,  is  used  in  simulating  the  tracking 
performance  of  the  Fire  Fly  gunsight  against  targets  in 
constant  altitude,  constant  airspeed,  constant  rate  turns 
at  bank  angles  of  30,  60,  and  70  degrees.  The  tracking  error 
responses  are  found  to  have  characteristics  similar  to  those 
observed  in  the  tracking  responses  obtained  from  man-in- 
the-loop  simulation  studies^onducted  by  McDonnell  Aircraft 
Company,  "^he  dependence  of  the  tracking  error  responses  on 
the  pilot  model  parameter  values  is  demonstrated  by  simulation 
runs  in  which  several  pilot  model  parameter  values  are  varied. 
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AN  INVESTIGATION  OF  THE  TRACKING  PERFORMANCE 


OF  THE  FIRE  FLY  MANUAL  DIRECTOR  GUNSIGHT 
FOR  AIR-TO-AIR  GUNNERY 

I.  Introduction 

Background 

Since  World  War  I,  the  air-to-air  gun  has  been  one  of 
the  fighter  pilot's  primary  weapons  in  aerial  combat.  Present 
guns  are  capable  of  unleashing  a lethal  amount  of  firepower 
somewhere  out  in  front  of  the  aircraft.  The  question  of  how 
best  to  control  and  use  that  firepower  with  some  type  of  fire 
control  system  is  currently  the  subject  of  considerable 
research  and  development  effort.  The  development  of  gunfire 
technology  was  essentially  dormant  from  195^-1964  as  a result 
of  the  emphasis  on  the  development  of  air-to-air  missile 
technology.  However,  recent  experience  in  fighter/fighter 
engagements  worldwide  indicates  that  where  sophisticated 
air-to-air  missile  systems  prevent  disengagement  from  an 
encounter,  the  outcome  will  probably  be  decided  by  a gunnery 
duel  (Ref  li  Chap.  II,  p.  1). 

Air-to-air  combat  between  present  day  jet  aircraft  may 
be  initiated  at  any  altitude  between  sea  level  and  50,000  feet, 
and  at  airspeeds  between  I50  knots  and  Mach  2+  (Ref  2i  5)» 
Aerial  combat  may  also  occur  between  aircraft  of  diverse 
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capabilities  and  multi-plane  engagements  are  common.  These 
factors  further  complicate  the  already  difficult  problem 
of  effective  fire  control. 


The  role  of  the  pilot  in  the  operation  of  present  fire 
control  systems  is  crucial.  In  most  systems,  the  pilot  is 
responsible,  at  least  in  part,  for  precision  tracking  and 
selection  of  the  appropriate  weapon  (missile,  rocket,  gun) 
to  be  used.  In  addition,  the  pilot  must  decide  when  to  engage 
or  disengage  an  enemy.  He  must  monitor  the  performance  of 
his  own  aircraft  and  he  must  continuously  be  aware  of  the 
threat  environment.  These  responsibilities  produce  a criti- 
cally high  pilot  workload. 

This  thesis  is  concerned  with  the  evaluation  of  the 
manual  director  gunsight  portion  of  an  integrated  automatic 
flight  control  and  fire  control  system  which  is  specifically 
designed  to  alleviate  pilot  workload  in  the  air-to-air  combat 
environment.  The  particular  system,  called  Fire  Fly,  was 
developed  by  the  Aircraft  Equipment  Division  of  General 
Electric  Company.  The  Fire  Fly  system  is  being  tested  in  a 
joint  program  conducted  by  the  Air  Force  Flight  Dynamics 
Laboratory  and  the  Air  Force  Avionics  Laboratory,  using  the 
F-I06  as  the  test  bed  aircraft.  The  overall  goal  of  the 
Fire  Fly  system  is  to  provide  integrated  automatic  flight 
control  and  fire  control  by  allowing  the  pilot  to  select 
from  a range  of  control  options.  By  selecting  the  option  he 
desires,  the  pilot  can  allow  the  Fire  Fly  fire  control  system 
to  automatically  control  the  aircraft  flight  control  system 
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for  precision  tracking  during  aerial  combat.  The  available 
range  of  control  options  extends  frcm  full  manual  control  by 
the  pilot  to  full  automatic  control  by  the  Fire  Fly  system. 

The  primary  objectives  of  the  Fire  Fly  system  are  to  provide 
improved  precision  tracking  and  to  reduce  pilot  workload 
during  aerial  combat. 

To  provide  the  reader  with  some  basic  insights  into 
the  fire  control  problem  and  its  solution,  the  following 
paragraphs  briefly  discuss  the  basic  fire  control  problem  and 
the  two  general  types  of  gunsight  system  which  can  be  used 
to  aid  the  pilot  in  the  precision  tracking  task. 

Stated  concisely,  the  fire  control  problem  for  air-to-air 
gunnery  isi  fire  a bullet  from  an  attacking  aircraft  at 
a target  aircraft  so  as  to  score  a hit  on  that  target.  To 
score  a hit,  the  solution  to  the  fire  control  problem  must 
be  obtained  and  displayed  to  the  pilot  in  a meaningful  and 
easily  used  form.  The  manner  in  which  the  solution  is  obtained 
is  a function  of  the  type  of  gunsight  system  used  (Ref  3i 
Chap.  I,  p.  1). 

In  the  highly  dynamic  environment  of  air-to-air  combat, 
the  attacker  can  seldom  point  his  gun  directly  at  the  target 
to  achieve  a hit.  To  obtain  a hit,  the  attacker  must  displace 
the  gun  direction  from  the  target  by  a lead  angle.  This  lead 
angle  is  a function  of  the  range  between  the  attacker  and  the 
target,  the  angular  rate  of  the  line-of-sight  from  the  attacker 
to  the  target,  target  motion,  and  bullet  ballistics.  The 
purpose  of  the  gunsight  is  to  provide  a measure  of  aiming 
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error  so  appropriate  firing  or  corrective  control  can  be 
effected  (Ref  1). 

All  air-to-air  gunsight  systems  direct  the  attacker 
gunfire  by  directing  the  attacking  pilot  to  fly  his  aircraft, 
and  hence  the  gun,  to  a particular  position  with  respect  to 
the  target.  This  task  is  accomplished  by  the  gunsight  system, 
which  displays  a dot  of  light  or  pipper  on  a Head  Up  Display 
(HUD)  in  front  of  the  pilot.  The  pilot  then  superimposes 
the  pipper  on  the  target  by  maneuvering  his  aircraft 
(Ref  3i  Chap.  VI,  p.  1). 

There  are  two  basic  classes  of  gunsight  system  for  air- 
to-air  gunnery  I the  disturbed  reticle  gunsight  and  the 
director  gunsight.  The  two  classes  differ  in  the  manner  in 
which  the  angular  rate  of  the  line-of-sight  from  the  attacker 
to  the  target  is  measured.  In  the  disturbed  reticle  system, 
the  lead  angle  is  computed  using  own  aircraft  angular  rate 
as  the  primary  input.  If  the  pilot  is  tracking  the  target 
so  that  the  pipper  remains  on  the  target,  then  own  aircraft 
body  rate  is  equal  to  the  line-of-sight  rate.  Thus,  the 
pilot  must  keep  the  pipper  on  the  target,  which  can  be 
dynamically  difficult  because  the  effect  of  feeding  back  own 
aircraft  angle  rates  to  the  sight  is  destabilizing  (Ref  5i  13-1^). 

The  commonly  used  Lead  Computing  Optical  Sight  (LCOS)  is  an 
example  of  a disturbed  reticle  sight. 

The  director  gunsight  uses  line-of-sight  rate  as  its 
major  input.  An  angle  tracking  radar  or  an  electro-optical 
tracker  is  used,  in  conjunction  with  a Kalman  filter,  to 
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estimate  the  angular  rate  of  the  line-of-sight  and  the 
target  acceleration  based  on  past  history  of  target  trajectory 
(Ref  2 1 8-9).  Although  the  director  system  is  more  complex 
and  expensive  than  the  disturbed  reticle  system,  it  offers 
distinct  advantages  for  air-to-air  combat.  The  director 
gunsight  provides  a direct  measure  of  gun  error  for  the  pilot 
to  null.  Also,  since  the  lead  angle  computation  does  not 
depend  on  own  aircraft  motion,  except  for  disturbance  errors, 
the  pilot's  control  task  is  not  complicated  by  sight  dynamics. 
In  addition,  if  angular  rate  measurement  noise  is  sufficiently 
attenuated  by  filtering,  the  director  system  provides 
improved  steering  stability  in  the  tracking  task  (Ref  5»  i9). 

The  Fire  Fly  system  employs  a director  gunsight  system 
which  utilizes  an  Ascot  electro-optical  tracker  and  a 
Kalman  filter  algorithm  for  estimation  of  the  target  state. 

Purpose 

The  purpose  of  this  thesis  effort  was  to  examine  the 
tracking  performance  of  the  Fire  Fly  manual  director  gunsight 
during  the  terminal  phase  of  air-to-air  combat.  The  evalu- 
ation was  conducted  using  a digital  computer  simulation. 

The  simulation  introduced  two  secondary,  but  important, 
objectives.  First,  through  the  development  of  the  digital 
computer  model  and  establishment  of  appropriate  simulation 
techniques,  the  author  established  the  usefulness  of  the 
ECS  Transient  Analysis  Integrated  Computer  Program  (EASY)# 
which  was  developed  by  Boeing  Computer  Services,  Inc,,  for 


analyzing  dynamic,  nonlinear  weapon  systems.  Second,  the  ] 

author  adapted  the  analytic  multi-axis  pilot  model,  which  was  il 

ij 

developed  in  the  Terminal  Aerial  Weapon  Delivery  Simulation  j 

(TAWDS)  program,  for  use  in  the  F- 10 6/Fire  Fly  gunsight  ( 

simulation.  | 

Scope 

Although  the  Fire  Fly  system  is  designed  for  use  in  air- 
to-air  gunnery,  air-to-ground  gunnery,  and  dive  bombing,  this 
thesis  was  limited  to  an  investigation  of  the  air-to-air 
gunnery  mission.  Specifically,  the  Fire  Fly  gunsight  system 
was  evaluated  only  in  the  manual  mode.  That  is,  the  tracking 
errors  were  presented  to  the  pilot  model,  which,  in  turn, 
commanded  the  aircraft  flight  control  system.  The  evaluation 
assumed  continuous  lock-on  of  the  tracker  and  perfect  estim- 
ation of  the  target  state.  Initial  conditions  were  specified 
to  put  the  attacker  at  an  altitude  of  10,000  feet  and  a speed 
of  Mach  .8.  Target  maneuvers  consisted  of  constant  altitude, 
constant  airspeed,  constant  rate  turns.  The  initial  position 
of  the  attacker  with  respect  to  the  target  was  specified  so 
as  to  put  the  attacker  near  a gun  solution.  The  attacker  was 
assumed  to  have  sufficient  available  thrust  to  maintain  air- 
speed throughout  the  tracking  maneuver. 

The  digital  computer  model  of  the  F-106  aircraft  and 
applicable  aerodynamic  data  were  provided  by  the  Air  Force 
Flight  Dynamics  Laboratory.  The  digital  computer  pilot  model 
was  adapted  from  the  pilot  model  used  in  the  Terminal  Aerial 


Weapon  Delivery  Simulation  (TAWDS)  programs  (Ref  6i  65-II9). 
Overview 

In  the  next  chapter  the  Fire  Fly  system  is  discussed. 

The  performance  objectives  of  the  Fire  Fly  system  are  pre- 
sented from  the  standpoint  of  the  integrated  flight  control/ 
fire  control  philosophy.  The  weapon  delivery  and  gunnery 
equations  are  derived,  and  expressions  are  developed  for 
elevation  and  traverse  tracking  errors  which  are  to  be  dis- 
played to  the  pilot.  In  addition,  a brief  discussion  of  the 
Kalman  filter  is  provided  from  the  viewpoint  of  the 
information  that  the  Kalman  filter  provides  to  the  gunsight. 

Chapter  Three  discusses  the  F-I06  aircraft  and  presents 
the  six  degree  of  freedom  equations  of  motion  which  describe 
the  aircraft.  In  addition,  the  longitudinal  and  lateral- 
directional  flight  control  systems  are  presented  and  discussed. 
In  Chapter  Four,  the  pilot  model  is  discussed.  The  generic 
model,  which  was  developed  for  the  F-4  aircraft  in  the  TAWDS 
program,  is  presented,  and  the  components  of  the  model  are 
discussed.  The  specific  model  used  for  the  F-I06  is  presented, 
and  justification  for  the  particular  choice  of  parameter 
values  is  given. 

Chapter  Five  presents  the  results  of  the  digital  computer 
simulation.  The  tracking  performance  of  the  closed  loop  system 
is  evaluated  against  targets  in  straight  and  level  flight, 
and  in  constant  altitude,  constant  airspeed,  constant  rate 
turns  at  bank  angles  of  30*  60,  and  70  degrees.  Time  history 


plots  of  elevation  and  traverse  tracking  errors  are  presented. 

Plots  of  corresponding  longitudinal  and  lateral  pilot  stick 
forces  are  also  presented.  The  changes  in  tracking  performance 
which  result  from  changing  the  nominal  pilot  model  parameter 
values  are  also  discussed.  Chapter  Six  concludes  that  the 
Fire  Fly  manual  director  gunsight  system  is  capable  of 
tracking  the  targets  simulated  and  that  the  EASY  computer 
program  is  useful  for  modeling  and  analyzing  dynamic, 

closed  loop  weapon  systems.  In  addition,  the  TAWDS  pilot  j 

model  can  be  adapted  for  use  in  the  EASY  computer  simulation  i 

of  the  F-106/Fire  Fly  gunsight  system.  Evaluation  of  the 

Fire  Fly  gunsight  tracking  performance  against  dynamically 

maneuvering  targets  is  recommended  as  is  validation  of  the 

digital  pilot  model  by  comparison  with  F-106/Fire  Fly 

man- in- the- loop  simulation  information. 

Appendix  A presents  a listing  of  the  EASY  computer  model 
used  in  the  simulation.  The  flow  of  the  model  is  discussed 
as  are  specific  modeling  and  simulation  techniques.  Appendix 
B contains  the  specific  aerodynamic  data  used  in  the  simulation. 
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II,  The  Fire  Fly  System 


General 

The  general  Fire  Fly  system  is  discussed  in  this  chapter 
to  give  the  reader  an  insight  into  the  type  of  fire  control 
system  which  is  being  considered  for  implementation  in 
future  aircraft  weapon  systems.  The  Fire  Fly  system  in- 
corporates the  concept  of  integrated  flight  control/fire 
control  by  providing  automatic  aircraft  control  to  aid  the 
pilot  in  the  weapon  delivery  task.  The  system  was  designed 
to  be  used  for  air-to-air  gunnery,  air-to-ground  gunnery,  and 
conventional  dive  bombing.  The  pilot  has  the  capability 
to  select  a variable  amount  of  automatic  control  authority. 

In  the  manual  mode,  the  pilot  performs  all  aircraft  control 
functions.  By  appropriately  selecting  authority  limits  for 
the  Fire  Fly  system,  the  pilot  relinquishes  a specified 
degree  of  responsibility  for  precision  tracking  control  to 
the  automatic  system.  Within  the  assigned  authority  limits, 
the  Fire  Fly  computer  provides  appropriate  tracking  control 
commands  to  the  aircraft  control  augmentation  system, 

A general  block  diagram  of  the  Fire  Fly  system  is  shown 
in  Figure  1,  The  tracking  system  provides  target  information 
to  the  Fire  Fly  computer.  The  computer  uses  this  information 
to  estimate  the  target  state.  Aircraft  sensors  provide  infor- 
mation about  own  aircraft  state  which  the  computer  compares 
to  the  target  state  estimate  to  generate  the  target  relative 
state.  The  target  relative  state  information  is  used  in  the 
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Figure  1.  Fire  Fly  System  Block  Diagram  (Adapted  from  Ref  li  Chap,  III, 


r 

weapon  delivery  calculations  to  compute  the  correct  gun 
solution.  The  correct  angles  are  then  compared  to  the  actual 
gun  position  relative  to  the  target  to  generate  error  signals. 
These  error  signals  are  displayed  to  the  pilot  on  a HUD 
(Ref  li  Chap.  IIIi  p,  3^)*  In  the  manual  mode  of  operation, 
the  pilot  attempts  to  null  the  errors  by  maneuvering  the 
aircraft.  When  the  Fire  Fly  system  is  given  automatic 
control  authority,  the  error  signals  are  also  fed  to  the  Fire 
Fly  control  laws  which  generate  commands  to  drive  the  control 
augmentation  system.  The  Fire  Fly  control  law  commands 
replace  pilot  tracking  control  commands  up  to  the  authority 
limits  selected  by  the  pilot.  Figure  2 depicts  the  flow  of 
information  in  the  Fire  Fly  computer. 

Since  the  investigation  reported  herein  concerns  only 
the  manual  mode  of  operation,  the  pilot  must  maneuver  the 
aircraft  to  null  the  gunnery  errors.  Furthermore,  since 
continuous  tracker  lock-on  and  perfect  estimation  of  the 
target  state  and  own  aircraft  state  are  assumed,  target  rel- 
ative state  information  is  presented  directly  to  the  weapon 
delivery  equations.  This  information  is  subject  to  a time 
delay  which  represents  the  time  required  for  the  filter  to 
process  target  and  own  aircraft  state  information. 

In  the  case  where  the  target  is  dynamically  maneuvering, 
the  delay  for  filter  processing  will  depend  on  the  number  of 
iterations  required  for  the  filter  to  attain  good  estimates 
of  the  target  state.  However,  for  the  cases  considered  in 
this  thesis,  where  the  target  acceleration  is  constant 
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Figure  2.  Fire  Fly  Computer  Software  Flow  Diagram 
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throu^jhout  each  simulation  run,  the  time  delay  is  due  only 
to  the  actual  information  processing  time.  Since  the  Kalman 
filter  has  an  iteration  interval  of  I/30  second  (Ref  7i  45), 
the  effects  of  the  information  processing  delay  are  negligible, 
as  the  simulation  results  in  Chapter  Five  indicate. 

Figure  3 depicts  a simplified  block  diagram  of  the 
closed  loop  tracking  task.  The  next  section  develops  the 
weapon  delivery  and  gunnery  equations  and  derives  expressions 
for  the  elevation  and  traverse  tracking  errors  which  are 
displayed  to  the  pilot. 

Weapon  Delivery  Equations 

The  weapon  delivery  equations  use  target  relative  state 
information  and  measured  attacker  parameters  to  calculate 
I the  gunnery  errors.  Figure  4 depicts  the  air-to-air  gunnery 

situation. 

The  relationships  shown  in  Figure  4 are  strictly  valid 
only  if  the  target  acceleration  is  constant.  If  the  target 
acceleration  is  time-varying,  the  relationships  depicted  in 
Figure  4 represent  an  approximation  to  second  order.  For  the 
development  of  the  weapon  delivery  equations  in  this  chapter 
the  second  order  approximation  is  assumed  to  be  valid. 

Referring  to  Figure  4,  R^  represents  the  present  required 
range  vector  for  a hit,  as  measured  from  the  attacker  to  the 
target,  R^^  represents  the  range  vector  from  bullet  release 
point  to  bullet  impact  point.  7^  is  the  attacker  velocity 
vector  relative  to  the  airmass,  and  7^  is  the  gun  muzzle 
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Figure  3*  Simplified  Closed  Loop  Tracking  Task 


Figure  4,  The  Air-to-Air  Gunnery  Situation 


velocity  vector.  is  the  target  velocity  vector  relative 
to  the  airmass,  and  A^  is  the  target  inertial  acceleration 
vector.  is  the  bullet  time-of-flight  and  g is  the  gravity 

vector.  The  present  required  range  vector  from  the  attacker 
to  the  target  is  given  by 


V T 
t^f 


- iA^T 


t"f 


(1) 


The  range  vector  which  represents  the  position  of  the  bullet 
after  one  time-of-flight  is  given  by 

'*ba  ' LV(''a  ^ (2) 

where  is  the  magnitude  of 
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Substituting  Eq  (2)  into  Eq  (1)  yields 


The  target  velocity  vector  and  the  target  acceleration 
vector  are  given  by 

(^) 


+ S 


(5) 


where  is  the  target  relative  velocity  vector,  is  the 
target  relative  acceleration  vector,  and  A is  the  attacker 
acceleration  vector.  The  attacker  acceleration  vector  can 
be  expressed  as 

A = A^  + g (6) 

where  A^  is  the  attacker  specific  force  vector,  measured  by 
the  outputs  of  body-mounted  accelerometers.  Substituting 
Eq  (6)  into  Eq  (5)  yields 


^t  “ -^r  ^b  S 


Substituting  Eqs  (4)  and  (?)  into  Eq  (3)  yields  the  range 
vector  required  for  a hit  as 

K ' tV(''a  * CV<''a  * ''n.)  - 

- 4(X^  + I^)t/ 


(?) 


(8) 
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The  ballistics  relationship  between  bullet  time-of-flight 
and  future  range  is 

Tf  • D/(V^  * (9) 


where  is  the  ballistic  coefficient  and  p/p^  is  relative 
air  density  (Ref  1:  Chap.  Ill,  p.  55). 


Letting 


b o^'  ^o  a m 


and  substituting  in  Eq  (9)  yields 


Solving  Eq  (11)  for  in  terms  6f  T^  yields 

“f  ' (''a 


Letting 


Kj,  = 1/(1  + Cj,Tj,) 


and  substituting  Eq  (I3)  into  Eq  (12)  yields 


“f  = ''b<''a  ^ \»s: 


Substituting  Eq  (14)  into  Eq  (8)  yields  the  range  vector 
required  for  a hit  as 


K “ ‘^b'^A  ■ ‘ ^f^r  ’ ^r^'^f' 
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Figure  5»  Body- Fixed  Coordinate  Systera 

Eq  (15)  represents  the  physical  vector  relationships 
between  the  attacker  and  the  target  with  respect  to  inertial 
space.  To  facilitate  the  use  of  these  relationships  for 
displaying  information  to  the  pilot,  the  inertial  vector 
quantities  in  Eq  (I5)  are  coordinatized  in  terms  of  an  axis 
system  fixed  in  the  body  of  the  attacking  aircraft.  Figure  5 
depicts  a body-fixed  coordinate  system,  u^,  where  u is  out 
the  nose  of  the  aircraft,  v is  out  the  right  wing,  and  w is 
nominally  down.  The  angle  Ej^^  is  the  angle  between  the  u 
axis  and  the  gun  axis.  The  angle  a is  the  aerodynamic  angle 
of  attacks  that  is,  the  angle  between  the  aircraft  velocity 
vector  and  the  u axis. 

Separating  Eq  (I5)  into  scalar  equations,  expressed  in 
terms  of  the  components  of  the  body-fixed  axis  system,  and 
assuming  that  the  component  of  in  the  v direction  is  zero, 
and  that  the  component  of  V in  the  v direction  is  zero, 

3L 
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yields  (Ref  7i  45) 

''ru  ' -V-'SV‘>=<'^Lg)  * V 4’'f<*bu  * \u> 

+ TjK^C^V^cos(a)]  (16) 

'’rv'  -T/''rv'^  <*’’> 

Rrw  ” * ''r»  i^f<V  *rw> 

+ TfK(,C(,V3Sin(a)]  (18) 

where  R i R ..i  and  R_  , are  the  components  of  the  ran^e  vector 
ru  rv  rw 

required  to  obtain  a hit,  expressed  in  aircraft  body-fixed 
coordinates. 

The  error  signals  to  be  displayed  to  the  pilot  can  be 
defined  as  the  angular  difference,  in  radians,  between  the 
range  components  required  for  a hit  and  the  present  measured 
range  components.  Thus  (Ref  7«  45) 


LV  ■■  ^^rw  ■ ^tw^/^ 

(19) 

LW  “ "^^rv  ■ ^tv^/^ 

(20) 

where  e^^y  and  are  the  gunnery  errors  in  elevation  and 
traverse,  respectively,  R^^  and  R^^  are  components  of  the 
present  target  range,  and  R is  the  magnitude  of  the  present 
range  vector.  Positive  elevation  error,  e^yt  is  decreased  by 
positive  pitch  rate,  and  positive  traverse  error,  is 

decreased  by  positive  yaw  rate.  Positive  pitch  rate  is 


defined  as  positive  rotation  about  the  body-fixed  v axis,  as 
given  by  the  right  hand  rule.  Positive  yaw  rate  is  defined 
as  positive  rotation  about  the  body-fixed  w axis,  as  given 
by  the  right  hand  rule. 

A measure  of  bullet  time-of-flight  is  required  for 
computing  the  range  vector  components  required  for  a hit. 

A method  of  calculating  bullet  time-of-flight  is  developed 
in  the  next  section. 

Time- of- Flight  Calculation 


The  scalar  relationship  between  future  target  range  and 
bullet  time-of-flight  is  illustrated  in  Figure  6. 


Figure  6.  Time-of -Flight  Calculation  (from  Ref  li  Chap. Ill, p. 54) 
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From  Figure  6, 


cos(X  - y)  » QR  + V^TjCOs(A^)VCDj  + iA^T^^sin(A^)3  (21) 


where 


Y = [vy(v^  * v„)3ag 

and  X is  the  sight  lead  angle,  A^  is  the  angle-off  of  the 
target  heading,  and  a is  the  gun  angle  of  attack.  Also 

O 


from  Figure  6, 


V^cos(A^)  = V^cos(X  - ttg)  + R 


where  R is  the  range  rate.  Substituting  Eq  (23)  into 
Eq  (21)  yields 

cos(X  - Y)  « [R  +Lv^cos(X  - ttg)  + R]T^}/[D^  + f A^T^^sin(A^) ] 
Solving  Eq  (24)  for  yields 

Df  = [R  + V^T^cos(X  - a ) + RTjVcos(X  - y)  - iA^T^^sin(AQ ) 


Recalling  Eq  (11),  which  is  repeated  here  for  convenience. 


Tf  ' V<''a  \ 


and  substituting  Eq  (25)  into  Eq  (11),  yields  the  following 
quadratic  equation  for  T^i 


C^[V^cos(X  - ttg)  + R - iA^T^^sin(A^)cos(X  - Y)3T^^ 

+ Cc,R  + R + V rcos(X  - a ) - cos(X  - y)I] 

u « ^ 

- ■*■  “ Y)3t^  + R = 0 (26) 

where  is  some  nominal  value  of  T^.  Since  the  effects  of 
target  acceleration,  angle-off,  lead  angle,  and  angle-of- 
attack  are  second  order,  nominal  values  are  used  in  the 
calculation.  Thus,  the  time-of-f light  computation  can  be 
mechanized  as 

+ R - 0^)t/  + iC^R  + R - C2)T^  + R = 0 (2?) 

where  Cj^  and  C2  are  constants  (Ref  Ij  Chap,  III,  p.  56). 

Eq  (27)  is  easily  solved  using  the  quadratic  formula  and  the 
smaller  value  is  chosen  as  the  bullet  time-of-flight. 

This  chapter  presented  a general  view  of  the  Fire  Fly 
integrated  flight  control/fire  control  system.  A simplified 
diagram  of  the  closed  loop  tracking  task  for  the  manual  mode 
of  operation  was  developed.  The  weapon  delivery  and  gunnery 
equations  were  developed  as  was  a method  for  calculating 
bullet  time-of-flight.  The  next  chapter  presents  the  six 
degree  of  freedom  equations  of  motion  which  describe  the 
F-IO6  aircraft.  The  equations  include  second  order  stability 
derivatives  and  the  effects  of  external  fuel  tanks.  Also 
presented  is  the  basic  F-IO6  flight  control  system. 
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IIIi  The  F-IO6  Aircraft 


I 

I 


The  F-I06  is  a single  engine  all-weather  interceptor 
built  by  General  Dynamic s/Convair.  Two  models  of  the  F-I06 
exist,  the  single  seat  F-I06A  and  the  two  seat  F-I06B. 

Although  designed  for  training  purposes,  the  F-I06B  can  also 
be  used  for  combat.  The  aircraft  modeled  in  this  simulation 
analysis  is  the  F-I06B  configured  with  two  v/ing-mounted 
external  fuel  tanks.  This  chapter  develops  the  general 
equations  of  motion  which  are  used  to  describe  the  aircraft 
and  discusses  the  F-I06  flight  control  system.  The  nonlinear 
aero-coefficients  were  obtained  from  Reference  8.  The  numer- 
ical values  of  all  stability  derivatives  are  presented  in 
tabular  form  in  Appendix  B.  These  tabular  values,  which  were 
obtained  by  curve  fitting  flight  test  data,  were  provided  by 
the  Air  Force  Flight  Dynamics  Laboratory.  The  reader  is  also 
referred  to  Appendix  A which  presents  and  discusses  the  Qxgital 
computer  model  used  in  the  simulation. 

Equations  of  Motion 

In  this  section,  the  equations  which  describe  aircraft 
motion  are  presented.  First,  coordinate  systems  and  sign 
conventions  are  discussed.  Second,  the  assumptions  upon  which 
the  equations  of  motion  are  based,  are  listed.  Third,  the 
equations  which  relate  the  orientation  of  the  aircraft  with 
respect  to  an  Earth-fixed  coordinate  system  are  presented. 
Finally,  the  equations  of  translational  motion  and  the 
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Fifcure  7,  Motion  of  the  Aircraft  with  Respect  to  Inertial 
Space  Resolved  into  Aircraft  Body-Fixed  Axes. 


equations  of  angular  motion  are  presented. 

Coordinate  Systems.  Figure  7 depicts  the  motion  of  an 
aircraft  with  respect  to  inertial  space.  Referring  to 
Figure  7»  is  the  total  linear  velocity  vector,  U,  V,  and 
W arc  the  components  of  the  total  velocity  vector  V^,  expressed 
in  terms  of  a body-fixed  coordinate  system.  The  body-fixed 
axis  system  u^  is  the  same  axis  system  that  was  illustrated 
in  Figure  5 of  Chapter  Two.  w is  the  total  angular  velocity 
of  the  aircraft  with  respect  to  inertial  space.  P,  Q,  and  R 
are  the  components  of  expressed  in  terms  of  the  body-fixed 
coordinate  system. 

Figure  8 shows  the  aerodynamic  forces  and  moments  which 

act  on  an  aircraft  in  flightt  again  expressed  with  respect  to 

the  body-fixed  coordinate  system.  F , F , and  F„  are  the 

X y z 
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Figure  8.  Aerodynamic  Forces  and  Moments  Acting  on  an 
Aircraft  in  Flight. 

aerodynamic  force  components,  and  L,  M,  and  N are  the  aero- 
dynamic moment  components. 

Assumptions.  The  development  of  the  aircraft  equations 
of  motion  is  based  on  the  following  assumptions i 

1.  The  aircraft  is  a rigid  body. 

2.  The  mass  and  mass  distribution  of  the  aircraft  are 
constant  over  the  time  interval  of  interest. 

3.  The  Earth  is  assumed  to  be  non-rotating  with  respect 
to  inertial  space  since  the  rotation  rate  of  the  Earth  is 
sufficiently  slow  so  as  to  be  negligible  in  the  time  period 
of  interest.  The  Earth  is  also  assumed  to  be  flat  since  the 
effects  of  the  Earth's  curvature  are  negligible  for  the  problem 
being  considered, 

4.  The  uw  plane  is  a plane  of  symmetry. 
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Picture  9.  Orientation  of  the  Body- Fixed  Coordinate 
System  with  Respect  to  an  Earth-Fixed 
Coordinate  System  (from  Ref  9i  Chap.  II,  p.  I5) 

Euler  Angles.  The  orientation  of  the  aircraft  body- 

fixed  coordinate  system  with  respect  to  an  Earth-fixed,  and 

thus  inertial  by  assumption  3,  coordinate  system  is  specified 

by  the  sequential  rotation  through  ^ (the  heading  angle),  0 

(the  pitch  angle),  and  $ (the  roll  angle).  Figure  9 depicts 

the  rotations  with  the  positive  directions  as  shown.  In 

Figure  9#  the  Earth-fixed  coordinate  system  X„Y„Z„  is 

assumed  to  have  its  origin  at  the  center  of  mass  of  the 

aircraft.  First,  the  X„Y_2„  system  is  rotated  about  the  2^ 

£ £ £ £ 

axis,  in  the  positive  direction  as  determined  by  the  right 
hand  rule,  through  the  angle  >?,  and  thus  becomes  Xj^Yj^Zj^. 
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Next,  the  system  is  rotated  about  the  axis  through 

the  ansrle  0,  yielding  the  coordinate  system  Finally, 

the  system  is  rotated  about  the  axis  through  the 

angle  The  resulting  coordinate  system,  uw,  is  the  body- 

fixed  coordinate  system  (Ref  9i  Chap.  II,  p.  16).  The 

• • • 

components  0,  and  $ shown  in  Figure  9 can  be  projected 
along  the  aircraft  body  axes  u,  v,  and  w to  obtain 


P = $ - 'P  sin(0) 

(28) 

• 

Q = 0 cos($)  + ? cos(0)sin($) 

(29) 

• • 

R = -0  sin($)  + ^ cos(©)cos(^>) 

(30) 

• • 

Eqs  (28),  (29),  and  (30)  can  be  solved  for  ’P,  0, 

obtain  (Ref  lOj  I3) 

• 

and  $ to 

• 

'P  = Q sin('I>)/cos(0)  + R cos($)/cos(0) 

(31) 

• 

0 = Q cos($)'  - R sin($) 

(32) 

$ = P + Q sin($)tan(0)  + R cos($)tan(0)  (33) 

Equations  of  Translational  Motion.  The  equations  of 
translational  motion  can  be  developed  utilizing  Newton’s 
Second  Law  which  states  that  the  time  rate  of  change  of 
linear  momentum  is  equal  to  the  sum  of  the  externally  applied 
forces.  The  scalar  equations  of  translational  motion  expressed 
in  terms  of  the  body-fixed  coordinate  system  are  (Ref  10 i 8) 
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X = m(u  + iciW  - RV) 


(3^) 


Y = m(V  + RU  - PW)  p5) 

Z = in(W  + PV  - QU)  (36) 

• • • 

where  Tn  is  the  mass  of  the  aircraft,  U,  V,  and  W represent 
the  time  rate  of  change  of  the  linear  velocity  components, 
and  X,  Y,  and  Z represent  the  sum  of  all  externally  applied 
forces  acting  along  the  body-fixed  axes  u,  v,  and  w,  respec- 
tively. 

The  summed  force  terms  X,  Y,  and  Z are  composed  of  the 
effects  of  gravity,  the  aerodynamic  forces,  and  the  thrust 
forces.  Thus, 


^ = ^x  ^^x  ^ Fx^ 

(37) 

Y - e,  + F + F 

y y-r 

(38) 

2 = g + F + F 

(39) 

where  g^,  g,  . and  g_  are  the  components  of  the 
* X y ' z 

gravity  vector 

F.,*  the  aerodynamic  force  components,  and 

X y z 

F„  , F,  • and  F are  the  thrust  components. 

*ip  Yijt  *ip 

The  gravity 

components  are  (Ref  10 1 I3) 

gx  - sin(©) 

(4o; 

g = mg  cos(0)sin($) 

(41) 

g_  = mg  cos(0)cos  ($) 

Z 

(42) 

where  g = 32.2  ft/sec^. 


* 
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The  thrust  side  force,  F 


is  assumed  to  be  zero  and 


the  thrust  components  F and  F can  be  expressed  as 


Fx^  = T cos(€^) 


Fz^  = T sin(€^) 


where  T is  the  magnitude  of  the  thrust  and  is  the  angle 

between  the  u body  axis  and  the  direction  of  the  thrust  vector. 
For  the  F-IO6,  = 0,  thus 


F„  = 0 (4( 

The  longitudinal  aerodynamic  force  terms  F and  F can 

A y 

be  expressed  in  terms  of  the  coefficients  of  lift  and  drag. 
Thus  (Ref  11:  243) 


Fj^  «=  sin(a)  - cos(a)j 


Fj,  •=  - qS[Cj^  cos(a)  + Cq  sin(a)] 


where  C and  C ^ are  the  dimensionless  coefficients  of  lift 
and  drag,  respectively.  The  aerodynamic  angle-of-attack 


a is  given  by 


a = arctan(W/U) 


S is  the  wing  reference  area  and  q is  the  dynamic  pressure. 
The  coefficient  of  lift  for  the  F-I06  can  be  expressed  as 
(Ref  81  34) 


29 


(50) 


C,  CL 


+ C, 


+ Cx  a + C, 


where 


is  the  lift  coefficient  for  a = 0 

is  the  lift  curve  slope 

is  the  tank  lift  coefficient  at  a = 0 
T 


the  tank  lift  curve  slope 


’E 


is  the  change  in  lift  coefficient  due  to  elevator 
E 

deflection 

is  the  elevator  deflection,  in  radians 


a is  the  angle-of-attack,  in  radians. 

The  drag  coefficient,  Cj^,  is  determined  as  a function  of 
Mach  number,  lift  coefficient,  and  elevator  deflection.  One 
method  of  determining  the  total  drag  coefficient,  and  the 
metliod  used  in  constructing  the  digital  model  of  the  F-106 
(see  Appendix  A),  is  a linear  interpolation  between  two  Mach 
numbers.  Given  two  curves,  C^,  and  Cp,  , which  represent 

JJt  Uo 


values  of  the  drag  coefficient  at  given  Mach  numbers,  each 
plotted  as  a function  of  lift  coefficient  and  elevator 
deflection,  the  total  drag  coefficient,  can  be  calculated 
as 


C 


D 


AM^ 


AMj^ 


(51) 
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where  is  the  difference  between  a hieh  reference  Mach 

number  and  the  actual  Mach  number,  for  example,  AM  =( .9-Mach) , 

n 

and  AMj^  is  the  difference  between  the  actual  Mach  number  and 

a low  reference  Mach  number,  for  example,  AM^=(Mach- . 8) , 

is  a constant  whose  value  depends  on  the  reference  Mach 

numbers,  and  and  Cp.  are  found  by  interpolation  for  a 
^2  ^1 

given  lift  coefficient  and  elevator  deflection. 

The  aerodynamic  side  force,  F , can  be  expressed  as 

y 

^ - \ln  * 


F is  given  by  (Ref  8i  3^) 
^Lin 


qSCC  3 (C  P + G R)b/2V„  + 6 C 6 ] (53) 

jg  jp  yp  y5  « ^6  “ 


where  P is  the  roll  angular  velocity,  in  radians/second 
R is  the  yaw  angular  velocity,  in  radians/second 
b is  the  reference  wingspan 

is  the  magnitude  of  the  total  velocity 
is  the  rudder  deflection,  in  radians 
6^  is  the  aileron  deflection,  in  radians 
C is  the  change  in  side  force  coefficient  due  to  3 

yp 

C is  the  change  in  side  force  coefficient  due  to  roll 

yp 

rate 


C„  is  the  change  in  side  force  coefficient  due  to  yaw 

yR 


rate 


C is  the  change  in  side  force  coefficient  due  to  a 
change  in  rudder  deflection 
C is  the  change  in  side  force  coefficient  due  to  a 
change  in  aileron  deflection. 


The  sideslip  angle,  g,  is  given  by 


P = arcsin(V/V^) 


F , the  nonlinear  component  of  the  side  force,  can  be 
^NL 

expressed  as  (Ref  8j  34) 


(5^) 


^NL 


= qS[(C  a + C )a 


yp  2 


a 


3 |3l 


+ qSlC  3 + C (BL)  P/2V  ] (55) 


where  |B|  is  the  absolute  value  of  the  sideslip  angle 

(BL)^  is  the  tank  buttock  line;  the  distance  measured 
in  the  uw  plane,  between  the  aircraft  water  line 
and  the  centerline  of  the  external  fuel  tank 


C is  the  change  in  C due  to  a change  in  a 

a ^ 

C.,  is  a second  order  effect  of  a on  C 
C..  is  the  change  in  C due  to  a change  in 


^3|3l 


'3 


C is  the  effect  of  the  external  fuel  tanks  on  C 

Yfi  i 

P rn 


is  the  effect  of  the  external  fuel  tanks  on  C 
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The  summed  force  terms,  X,  Y,  and  Z can  therefore  be 


expressed  as 


X = -mg  sin{©)  + qSLCj^  sin(a)  - Cj^  cos(a)3  + T (56) 


Y = mg  cos(0)sin(^>)  + q3[C^  g + (C^  P K)b/2V 

Yp  Yp  Yp  A 


C 6p  + C 6.]  + qS[(C  a + C )a 


+ . lgl3P  + qSLC^  g -»■  (BL)„  P/2V^3  (57) 

^P|P|  ""t 


Z = mg  cos(0)cos($)  - qSLCT  cos(a)  + sin(a)I!  (58) 


Equations  of  Angular  ’'■"otion.  Like  the  equations  of 
translational  motion,  the  equations  of  angular  motion  are 
developed  by  applying  Newton's  Second  Law  to  angular  momentum. 
Newton's  Law  states  that  the  time  rate  of  change  of  angular 
momentum  is  equal  to  the  sum  of  the  externally  applied 
moments.  Thus,  the  scalar  equations  of  angular  motion, 
expressed  in  terms  of  the  body-fixed  coordinate  system, 
are  (Ref  lOi  11) 

L = PI^  + QR(I^  - ly)  -(PQ  + R)I^^  (59) 

M = Qly  + PR(I^  - I^)  (P^  - 

N = RIj^  + PQdy  - I^)  + (QR  - P)I^^  (61) 
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where  I^i  I„i  and  I are  the  moments  of  inertia  about  the  u, 

X y z 

V,  and  w axes 

is  the  product  of  inertia  in  the  uw  plane 

• t • 

P,  Q,  and  R are  the  time  rates  of  change  of  the  angular 
velocity  components 
L is  the  total  rolling  moment 
M is  the  total  pitching  moment 
N is  the  total  yawing  moment. 

The  total  rolling  moment i L.  can  be  expressed  as  the  sum 
of  a linear  component,  ^ nonlinear  component, 

where  the  latter  contains  the  effects  of  the  external  fuel 
tanks,  can  be  written  as  (Ref  8i  34) 


ll,in  = 53b[C  0 + (C  P C R)b/2V^  + C + C 6^3 


R 


R 


^z^^ref  - 


eg' 


y'^'ref  “ ^cg^ 


(62) 


where  b is  the  reference  wingspan 

F_  is  the  aerodynamic  force  in  the  w direction 
F is  the  aerodynamic  side  force 

y 

- Y is  the  eg  moment  arm  in  the  v direction 

ref  eg 

is  the  eg  moment  arm  in  the  w direction 

is  the  change  in  for  a change  in  sideslip 

3 

C.  is  the  change  in  C»  for  a change  in  roll  rate 
P * 

Cf  is  the  change  in  C»  for  a change  in  yaw  rate 
R ^ 
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is  the  change  in  for  a change  in  rudder 

6, 


'R 


deflection 


is  the  change  in  for  a change  in  aileron 
^ A 

deflection 

The  nonlinear  rolling  moment,  is  given  by 

(Ref  8i  3^-35) 


^NL  “ 5SbCpJ[(C»  Co  )a  + Co  I ] + [(C«  a 


u.  ■ v->A  < Uo  IK  I J ' L \ 0 

^0^2  |g|' 


+ Cjg  )P  + R>(b/2V^)}  + qSb[C^  pJ  - [C^  (BL) 


'P„  ""R 

a a 


+ C (2.33  ft)](BL)„  P/2V  J 


where  is  the  change  in  for  a change  in  a 


Bt  0|P 


(63) 


C,  is  a second  order  effect  of  a on  C 


-e.. 


C,  is  the  change  in  C»  due  to  a change  in  1(31 

^3  131  S - 1 I 

C.  is  the  change  in  C»  due  to  a chance  in  a 

Pa  -“P 

C,  is  a second  order  effect  of  a on  C- 

PaP  "P 

C.  is  the  change  in  C,  due  to  a change  in  a 

\ 'r 

C,  is  the  effect  of  the  tanks  on  C. 

^ rt  X 


C,  is  the  effect  of  the  tanks  on  C« 
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C is  the  effect  of  the  tanks  of  C 

The  total  pitching  mo’nent,  M,  is  the  sum  of  the  pitching 
moment  due  to  thrust,  a linear  pitching  moment,  "^j^in’ 


and  a pitching  moment  due  to  the  external  fuel  tanks,  M 
The  pitching  moment  due  to  thrust  is  (Ref  12 i 2?) 


Mm  = Z F - X F 
T o x^  o 


Tank* 


(64) 


but  since  F_  =0  from  Eq  (46) 
^T 


M-  = T 


(65) 


► 

f 


where  is  the  distance,  measured  in  the  w direction, 
between  the  aircraft  center  of  gravity  and  the  origin  of 
the  thrust  vector. 

The  linear  pitching  moment,  is  (Ref  8i  35) 


M 


Lin 


qSc[C^  + (Q  + a)c/Vm 
0 ‘ q 

- F2,[(NP)c  - + F^(Z 


+ C, 


m 


ref 


^cg) 


(66) 


where  c is  the  mean  aerodynamic  chord 

a is  the  time  rate  of  change  of  anerle-of-attack 
6g  is  the  elevator  deflection,  in  radians 
(NP)  is  the  aircraft  neutral  point 

X is  the  distance,  measured  in  the  u direction,  from 

the  wing  leading  edge  at  the  wing  root  to  the  aircraft 
center  of  gravity 


1 

i 
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L 


is  the  pitching  moment  coefficient  for  a ^ 0 
”o 

is  the  pitching  moment  coefficient  due  to  pitch 

q 

rate 

C is  the  pitching  moment  coefficient  due  to  elevator 
deflection. 

The  tank  pitching  moment  is  siiven  by  (Ref  8i  35) 


'"Tank  ' 5so(0^  ♦ a)  (6?) 

0 ^ 


where  and  are  the  pitching  moment  coefficients  due 

djp  oc^ 

to  the  external  tanks. 

The  total  yawing  moment,  N,  is  expressed  as  the  sum  of 
a linear  yawing  moment,  and  a nonlinear  yawing  moment, 

^NL’  includes  the  effects  of  the  external  tanks. 

The  linear  yawing  moment  is  given  by  (Ref  8t  35) 

"Lin  = 521>[C  B * (C  P * C R)b/2Vj  * °n 

B P R “r  ‘a 


+ F r(NP)c  - X 3 - - Y ) 


(68) 


where  C is  the  change  in  yawing  moment  coefficient  due  to 

"fj 

C is  the  change  in  C due  to  roll  rate 
np  n 

C is  the  change  in  C due  to  yaw  rate 
np  n 


is  the  change  in  due  to  rudder  deflection 
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C is  the  change  in  C due  to  aileron  deflection. 

n 


The  nonlinear  yawing;  moment  is  expressed  as  (Ref  8i  35) 


" qSb[(C„  a + C„  )a  + (C 


n |3|  )fi 

"did  I "d 


+ (C„  a + C, 


'-6  a"  ^-6  6 ^E)^:i-5Sb[C  3 

a“  °A®E  3™ 


- 0.6(BL)^  P/2V^] 


(69) 


where  C is  a second  order  effect  of  a on  C 
1 1 


n 


*3  2 

C is  the  change  in  C due  to  a 
"3c  "d 

^-3|3(  change  in  due  to  jfjj 

AC  is  the  incremental  chanere  in  C due  to  3 
“3  n 

is  the  change  in  C„  due  to  a 


n, 


n 


6 6 
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is  the  change  in  C due  to  6 


n, 


The  total  ’noments,  Lt  M,  and  N can  then  be  expressed  as 


L = qSbCc,  3 + (C,  P + C,  R)b/2V_  ^ C,  6 + C,  6 ] 

Xg  Xp  Xp  T X^^  R X^^  A 


^ ^z^^ref  - - ^y<^ef  ' ^cg)  ^ qSbC3[(C  a + C )a 


R]a(b/2V,p)} 


3 13  I 
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+ qSbCc,  3 - [C„  (BL)_  + (2.33  ft)](3L)T,  P/2Vj 

3^  0^  ^ % T T 


(70) 


M = qSc[C^  + (Q  + a)c/V^  + 6^]  - F^[(NP)c  - X ] 


0 q 


^ Px^^ef  - 2cg)  ^ ^ T + q3c(C^  + a)  (71) 


N = qSb[C„  3 + (C„  P + C„  R)b/2V^  + 6p  + 6^] 

^R 


np-  ' -np-  np 


* PyC(NP)o  - q3bL(C  c + 0 )a 


+ (C_  |3l  + AC  )3  + (C„  a + C„  6„)6.3 
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+ qSbCC^  3 - 0.6(BL)j  P/2V^3 
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(72) 
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In  this  section  the  six  de^^ree  of  freedon  nonlinear 
equations  of  aircraft  motion  were  developed.  These  equations 
account  for  the  effects  of  the  aerodynamic  forces  and  moments 


acting:  on  the  aircraft,  the  force  of  gravity,  the  effects 
of  thrust,  and  the  effects  of  external  fuel  tanks.  The 
equations  which  specify  the  orientation  of  the  aircraft  with 
respect  to  an  Earth-fixed  coordinate  system  were  also  presented. 
The  development  of  the  equations  of  motion  was  based  on  the 
assumptions  that  the  aircraft  is  a rigid  body  and  that  the 
aircraft  mass  and  mass  distribution  remain  constant  over  the 
time  interval  of  interest.  Further,  it  was  assumed  that  the 
Earth-fixed  coordinate  system  is  nonrotating  with  respect  to 
inertial  space,  and  finally,  that  the  body-fixed  uw  plane 
is  a plane  of  symmetry.  The  equations  describing  the 
orientation  of  the  aircraft  with  respect  to  the  Earth-fixed 
system  are  given  by  Eqs  (31)i  (32)»  and  (33)-  equations 

of  translational  motion  are  given  by  Eqs  (3^)»  (35) i and  (36), 
where  the  summed  force  terms  X,  Y,  and  Z are  given  by  Eqs 
(56),  (57)  • Sind  (58)1  respectively.  The  equations  of  angular 
motion  are  given  by  Eqs  (59).  (60),  and  (6I),  where  the  total 
moments  L,  M,  and  N are  given  by  Eqs  (70),  (71),  and  (72), 
respectively. 

The  above  equations  were  used  to  describe  the  motion  of 
the  F-I06  aircraft  in  the  digital  computer  simulation.  In 
order  to  completely  describe  the  aircraft,  one  must  specify 
not  only  the  equations  of  motion,  but  also  the  aircraft 
flight  control  system. 
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F~106  Flight  Control  System 


The  F-106  utilizes  elevens  for  both  pitch  and  roll 
controlt  and  a conventional  rudder  for  directional  control. 
The  elevens  deflect  symmetrically  for  pitch  attitude  changes 
and  differentially  to  initiate  roll  angle.  The  longitudinal 
flight  control  system  utilizes  a pitch  damper  as  shown  in 
Figure  10.  Pitch  rate  q is  fed  back  through  a scheduled 


gain  Kq  and  a washout  circuit  with  a scheduled  time  constant 

T . The  scheduled  gain  and  tim.e  constant  are  functions  of 

q 

the  air  density  ratio.  The  specific  values  are  listed  in 
Appendix  B.  The  output  of  the  washout  circuit  is  fed  through 
an  authority  limiter  with  limits  of  + 1 degree,  .and  then  to 
a series  actuator.  The  resulting  output  is  compared  to  the 


Figure  10.  F-106  Longitudinal  Flight  Control  System  (from 
Ref  13i  51). 
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pilot's  pitch  command,  6gp.  The  differenced  quantity,  6^^, 

is  fed  through  a pitch  limiter  which  limits  the  commanded 

input  to  between  +8  degrees  and  -25  degrees.  The  output  of 

the  limiter  is  fed  to  the  power  actuator  and  the  resulting 

output,  6g,  is  the  elevon  deflection  for  pitch,  in  degrees. 

The  lateral-directional  control  system  includes  a yaw 

damper  and  differential  elevon-to-rudder  crossfeed  for  turn 

coordination.  A block  diagram  of  the  lateral-directional 

control  system  is  shown  in  Figure  11.  The  yaw  damper  employs 

feedback  of  roll  rate  P and  yaw  rate  R.  The  roll  rate  is  fed 

through  a scheduled  gain  Kp  and  then  through  a simple  washout 

circuit.  The  scheduled  gain,  Kp,  is  a function  of  dynamic 

pressure.  Specific  values  of  Kp  are  listed  in  Appendix  B. 

The  yaw  rate  R is  fed  through  a gain  of  -1  and  then  a simple 

washout  circuit.  The  pilot's  differential  elevon  command, 

6^p,  is  fed  through  a power  actuator  and  then  to  the  elevon 

control  surfaces  as  6^.  The  commanded  roll,  is  also 

fed  through  a gain  K,  , and  added  to  the  pilot's  rudder 

® A 

command,  6pp,  The  gain  K^  is  scheduled  as  a function  of 

A 

dynamic  pressure.  The  specific  values  of  K^  are  also  listed 

°A 


in  Appendix  B. 

This  chapter  presented  the  equations  of  motion  which 
describe  the  F-106b  aircraft  configured  with  external  fuel 
tanks.  The  basic  flight  control  systems  were  also  presented. 
The  next  chapter  develops  the  analytical  pilot  model  which 
was  used  for  the  closed  loop  simulation. 


\ 

\ 
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Figure  11.  F-106  Lateral-Directional  Control  System  (from  Ref  I3i5l). 


IV.  The  Pilot  Model 

General  Considerations 

For  closed  loop  digital  simulation  of  the  air-to-air 
ffunnery  tracking  task,  a suitable  analytical  model  of  the 
pilot  must  be  developed.  This  analytical  model  should  relate 
the  air-to-air  gunnery  errors  to  aircraft  steering  commands 
in  a way  which  approximates  the  response  of  a human  pilot. 

To  realistically  develop  such  a model,  one  needs  to  know 
how  the  human  pilot  performs  in  relation  to  the  particular 
aircraft/weapon  delivery  system  and  the  specific  aerial 
gunnery  task  being  considered. 

Perhaps  the  most  suitable  method  for  obtaining  a measure 
of  the  human  pilot's  response  is  man-in-the-loop  simulation. 
From  man-in-the-loop  studies,  one  can  obtain  time  histories 
of  the  tracking  error  response  and  a measure  of  the  pilot's 
frequency  response  characteristics.  By  controlling  the 
tracking  task  so  that  it  can  be  duplicated  in  an  all  digital 
simulation,  one  seeks  to  develop  an  analytical  model  which 
exhibits  tracking  and  frequency  response  characteristics 
which  are  similar  to  those  observed  for  the  human  pilot.  For 
the  system  simulated  in  this  report,  however,  appropriate 
man-in-the-loop  studies  have  not  yet  been  completed.  There- 
fore, an  appropriate  analytical  pilot  model  must  be  developed 
from  other  sources. 

In  the  development  of  the  Terminal  Aerial  Weapon 
Delivery  Simulation  (TAWD3)  programs,  McDonnell  Aircraft 
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Company  (McAir)  proposed  and  validated  a multi-axis 
analytical  pilot  model  for  the  F-4  aircraft  for  use  in  aerial 
ffunnery  tasks,  A complete  discussion  of  the  TAWDS  pilot 
model  development  can  be  found  in  Reference  14.  The  TAWDS 
pilot  model  was  used  successfully  in  the  simulation  of  F-4 
aircraft  configurations  with  widely  varying  flying  qualities. 
Therefore,  it  was  felt  that  the  basic  structure  of  the  TAWDS 
pilot  model  was  appropriate  for  preliminary  evaluation 
of  the  F-106/Fire  Fly  director  gunsight  system. 

The  Generic  Pilot  ^'!odel 

The  TAWDS  pilot  model  for  air-to-air  gunnery  was  developed 
for  use  in  the  terminal  phase  of  air-to-air  tracking,  where 
the  tracking  task  was  assumed  to  be  deterministic  in  that 
the  target  was  not  maneuvering  in  a random  manner.  The 
result  of  this  assumption  is  that  the  pilot  directly  observes 
the  effect  of  his  control  inputs  on  the  tracking  errors 
(Ref  6i  66).  McAir  found  that  the  characteristics  of  the 
tracking  error  time  histories  for  man-in- the-loop  simulations 
were  similar,  regardless  of  the  particular  pilot,  the  aircraft 
flying  qualities,  or  the  characteristics  of  the  gunsight. 

The  elevation  tracking  error  contains  two  predominant 
frequency  components,  both  of  which  exhibit  a limit  cycle 
type  of  response.  The  high  frequency  component  is  due  to 
the  pilot's  interaction  with  the  aircraft  short  period 
dynamics,  and  the  low  frequency  component  results  from  the 
pilot's  interaction  with  the  closed  loop  poles  which  result 
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from  the  relative  geometry  between  the  attacker  and  the  target. 
The  traverse  tracking  error  is  a single  frequency  limit  cycle 
response  with  a period  of  four  to  eight  seconds  (Ref  6t  68-69). 
The  TAWDS  pilot  model  was  structured  so  as  to  exhibit  responses 
in  the  elevation  and  traverse  tracking  error  time  histories 
which  were  similar  to  those  observed  for  the  human  pilot. 

A block  diagram  of  the  TAWDS  pilot  model  is  shown  in 
Figure  12.  The  model  represents  the  pilot  as  a proportional 
plus  derivative  observer  of  both  the  elevation  and  traverse 
tracking  errors,  with  deadzones  on  the  error  rates.  The 
simple  lag  networks  which  filter  the  true  errors  represent 
the  pilot's  response  delay  and  subsequent  smoothing  of  errors 
which  contain  noise  or  which  have  frequencies  above  the 
pilot's  observation  bandwidth  (Ref  6:  71). 

In  the  longitudinal  portion  of  the  pilot  model,  the 

smoothed  error  signal  e^,.^,  is  fed  through  a deadzone  DZ, 

LVS 

The  outputs  of  the  deadzone  and  the  proportional  gain  block 
Kp  are  summed  to  produce  a projected  error  This 

projected  error  represents  the  high  frequency  component  of 
the  tracking  error.  The  smoothed  error  signal  is  also 

fed  through  a second  order  low  pass  filter  and  a gain  Kj^. 

The  resulting  error  signal  is  the  low  frequency  component  of 
the  tracking  error.  This  low  frequency  component  is  then 
added  to  the  projected  error  resulting  signal  is 

multiplied  by  the  pilot's  transmission  gain,  io  determine 

the  pilot's  rate  input  to  the  control  stick,  6p  (Ref  6i  70-72). 
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Error  (Rad) 


Fiarure  12.  Generic  f'lulti-Axis  Pilot  Model  for  Air-to-Air  Gunnery  (from  Ref  14i  2?) 


In  the  lateral  portion  of  the  pilot  nodel,  the  smoothed 
error  is  fed  throusch  a proportional  sjain  block,  K__d. 

Lifio  KIK 

The  smoothed  error  siarnal  is  also  differentiated  and  fed 


throueih  a deadzone  DZ' 


LW 


The  output  of  the  deadzone  is 


co’Tibined  with  the  sieinal  fro’n  the  Kpjrjpj  block  to  form  the 
projected  lateral  error  This  signal  is  added  to  the 

signal  from  the  output  of  the  gain  block  K^.  The  portion 
of  the  model  which  feeds  differential  bank  angle  through  the 
deadzone  DZ^  and  the  gain  represents  the  pilot's  ability 
to  adjust  his  bank  angle  based  on  his  perception  of  the 
bank  angle  difference  between  his  aircraft  and  the  target. 
After  combining  the  projected  error  and  the  error  due  to 
differential  bank  angle,  the  resulting  signal  is  fed  through 
the  pilot  lateral  transmission  gain,  K, is  shown  as 


"AIL* 
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a function  of  elevation  tracking  error,  e^yi  since  the  pilot 
tends  to  act  as  a radial  controller  in  attempting  to  null 
both  elevation  and  traverse  errors  simultaneously.  The 
projected  error  is  also  fed  through  a rudder  pedal  gain 
since  the  pilot  will  probably  try  to  coordinate  his  lateral 
stick  and  rudder  commands  (Ref  6i  7^-76). 

The  TAWDS  studies  determined  a range  of  parameter  values 
for  the  multi-axis  pilot  model.  Table  I contains  a listing 
of  the  range  of  recommended  values  for  an  F-4  aircraft 
configured  with  a director  gunsight.  Because  the  parameters 
listed  in  Table  I affect  the  stability  of  the  closed  loop 
tracking  task,  the  specified  values  are  strictly  applicable 


1 
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only  for  the  F-4  aircraft.  However,  the  recommended  range 
of  values  was  used  as  a guide  in  choosing  first  trial 
parameter  values  for  the  F-IO6  pilot  model.  The  following 
paragraphs  discuss  the  function  of  each  parameter  in  the 
pilot  model.  This  discussion  is  adapted  from  Reference  14, 
pages  29  through  34. 


Table  Pilot  !'-'!odel  Parameter  Values  for  an  F-4  Aircraft 
with  Director  Gunsight  (from  Ref  FT  73,  7^) 


Longitudinal  Kpp  K„ 

Parameters  ^ ” 

Range  of  Values  0.75  ^CO  *75-1.25 


1. 0-3.0 


Longitudinal  t- 

Parameters  * 

Range  of  Values  O.O-O.05 


<%• 


1.0,  0.6 


5.0 

mr/sec 


Lateral 

Parameters 

Range  of  Values 


Lateral 

Parameters 

Range  of  Values 


0. 0-0.1  1.0 


<.36 


0.0-0.05  5.0 

mr/sec 


Each  parameter  in  the  pilot  model  contributes  to  the 
overall  objective  of  matching  as  closely  as  possible  the 


tracking  characteristics  of  the  human  pilot.  The  prefilter 
time  constants,  Tg  and  T,p,  represent  the  pilot's  ability  to 
smooth  observed  errors.  The  frequency  and  damping  parameters 
of  the  pilot's  low  pass  filtering  process,  and  repre- 
sent the  pilot's  interaction  with  the  relative  geometry  of 
the  elevation  tracking  task.  The  gain  determines  the 
amplitude  of  the  low  frequency  component  in  the  elevation 
tracking  error.  The  deadzone  on  the  elevation  error  rate, 

DZ*  , produces  the  desired  limit  cycle  response  in  the  high 
®LV 

frequency  component  of  elevation  tracking  error.  The  gain 
Kq  affects  the  amplitude  of  the  high  frequency  component. 

The  longitudinal  pilot  gain,  Kpg,  also  affects  the  amplitude 
of  the  high  frequency  component.  Referring  to  Table  I,  it 
can  be  seen  that  the  range  of  recommended  values  for  Kpg 
is  a function  of  is  the  value  of  gain  at  which 

the  closed  loop  system  poles  due  to  the  aircraft  short 
period  mode  pass  from  the  left  half  s-plane  to  the  right 
half  s-plane.  Thus  the  value  of  Kpg  chosen  must  be  large 
enough  to  excite  the  high  frequency  mode  of  the  elevation 
tracking  error,  but  not  so  high  as  to  cause  the  tracking  to 
become  unstable.  The  value  of  Kpg  is  also  a function  of 
the  target  range,  but  constant  gain  settings  were  used  in 
the  TAWDS  program  development.  In  the  F-IO6  simulation 
studies  described  in  this  report,  constant  values  were 
chosen  for  Kp^  since  all  simulation  runs  were  begun  at  the 
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same  range,  and  the  overall  change  in  range  magnitude  was  not 
significant  enough  to  warrant  gain  scheduling. 

The  deadzone  DZ*  describes  the  pilot's  ability  to 
®LW 

perceive  traverse  error  rates.  The  gain  parameters 
and  affect  the  amplitude  and  frequency  of  the  traverse 

error.  The  values  of  and  are  set  so  as  to  stabilize 

the  lateral  tracking  task  and  produce  the  desired  four  to 
eight  second  period  response.  shown  as  a function 

of  true  elevation  tracking  error  to  represent  the  pilot's 
intentional  coupling  as  he  tries  to  null  elevation  and  traverse 
errors  simultaneously.  As  with  Kpg,  constant  values  of 
were  used  in  both  the  TAWDS  and  F-IO6  simulation  studies. 

The  gain  param.eter  models  the  pilot's  rudder  pedal 

commands.  Man-in- the-loop  simulations  by  McAir  showed  that 
the  pilot's  rudder  force  commands  were  usually  negligible 
during  terminal  air-to-air  tracking.  The  deadzone  and 
the  gain  parameter  relate  pilot  stick  forces  due  to 
the  sensed  differential  bank  angle  between  the  target  and 
attacker.  Differential  bank  angle  feedback  is  probably  not 
a primary  input  during  terminal  gunnery  tracking,  but  the 
pilot  may  use  differential  bank  angle  cues  to  prevent 
excessive  wing  rock  (Ref  14 « 32). 

McAir  was  successful  in  utilizing  the  multi-axis 
analytical  pilot  model  discussed  above  to  obtain  tracking 
responses  similar  to  those  of  human  pilots  performing  the 
same  task.  The  TAWD3  studies  tested  the  pilot  model  with 
F-4  aircraft  models  which  were  configured  so  as  to  produce 
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widely  varyine:  sets  of  flying;  qualities.  Because  the  TAWDS 
pilot  model  was  used  successfully  reffardless  of  the  aircraft 
flying;  qualities,  the  same  basic  structure  was  chosen  for  the 
F-IO6  digital  simulations  performed  in  the  evaluation  of  the 
Fire  Fly  director  gunsight.  A discussion  of  the  specific 
F-IO6  pilot  model  is  presented  in  the  next  section. 


F-IO6  Pilot  Model 


The  structure  of  the  F-IO6  pilot  model  is  basically  the 

same  as  the  TAWDS  pilot  model,  and  parameter  values  were 

chosen  to  be  nominally  within  the  range  of  values  recommended 

in  Table  I.  The  value  of  Kp^,  the  longitudinal  pilot  gain, 

was  chosen  by  root  locus  analysis.  This  value  was  chosen  to 

be  high  enough  to  stabilize  the  low  frequency  component  of  the 

elevation  tracking  error.  The  upper  value  of  was  limited 

so  as  not  to  cause  the  complex  poles  which  are  due  to  the 

pilot's  interaction  with  the  aircraft  short  period  dynamics, 

to  cross  over  into  the  right  half  s-plane  (Ref  14:  55)-  If 

the  deadzone  DZ*  is  removed  from  the  pilot  model,  the 
®LV 

longitudinal  pilot  model  transfer  function  from  elevation 
tracking  error  to  longitudinal  force  command  is  (Ref  14:29) 


Kp.[s-^  + (2^u>  +Kp)s^  + (O'  ^+K^2<  Js  + (K,-*-Kp)u.'/] 


(73) 
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The  parameters  of  the  transfer  function  were  set  to 


nominal  values  and  a root  locus  for  the  closed  loop  longi- 
tudinal system  was  determined  as  a function  of  (The 

EASY  Analysis  Program  provides  a straightforward  tool  for 
producing  root  loci  about  a linearized  operating  point). 

A partial  root  locus  for  the  longitudinal  system  is  shown 
in  Figure  I3.  The  longitudinal  system  has  two  open  loop 
poles  which  start  at  s = 0 and  move  into  the  right  half  s-plane 
and  then  back  toward  the  pilot  model  complex  zeros  as  the 
gain  Kpg  is  increased.  For  the  nominal  pilot  model  parameters 


Figure  I3.  Partial  Root  Locus  for  Closed  Loop  Longitudinal 
System. 
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Figure  14.  Partial  Root  Locus  for  Closed  Loop 
Lateral  System. 

the  minimum  value  of  Kp^  for  stability  is  Kpg  = 0,5,  As 
Kpg  is  increased,  the  complex  poles  due  to  the  aircraft 
short  period  mode  move  toward  the  imaginary  axis.  The  value 
of  Kpg  for  which  the  root  locus  branch  due  to  the  aircraft 
short  period  crosses  the  imaginary  axis  is  the  crossover 
gain,  From  Figure  I3,  this  value  is  Kp^  = » 1,5. 

The  nominal  gain  value  chosen  for  the  simulation  analysis  was 


In  the  lateral  pilot  model,  the  value  of  the 

AIL 

lateral  pilot  transmission  gain,  must  be  set  high  enough  to 

stabilize  the  complex  poles  resulting  from  the  relative 

geometry  of  the  encounter  (Ref  14 1 55).  By  removing  tlie 

deadzone  DZ«  from  the  lateral  pilot  model  and  assuming  K 
®LW  «F 
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and  Kpyjj  to  be  zero#  the  transfer  function  from  traverse 
tracking  error  to  lateral  force  command  is 


^LAT  ^RTR^ 


(74) 


'LW 


trpS  1 


with  the  parameters  set  to  nominal  values,  a root  locus 
for  the  closed  loop  lateral  system  was  determined  as  a 
function  of  A partial  root  locus  for  the  lateral 

system  is  shown  in  Figure  14.  The  open  loop  rela'';ive 
geometry  poles  at  s = 0 move  into  the  right  half  s-plane  and 
then  back  toward  the  zero  at  s = as  the  value  of 

is  increased.  The  minimum  value  of  for  system  stability 

is  = 0.2.  For  the  simulation  analysis,  a value  of 

= 0.36,  which  corresponds  to  the  maximum  recommended 
value  in  Table  I,  was  used  as  the  nominal  value. 

Figure  I5  is  a block  diagram  of  the  pilot  model  used 
in  the  F-106/Fire  Fly  director  gunsight  simulation.  The 
parameter  values  shown  in  each  block  of  the  figure  are  the 
nominal  values  which  were  used  for  the  simulation.  Two 
adjustments  were  made  to  the  pilot  model  structure.  First, 
the  contribution  of  the  differential  bank  angle  portion 
of  the  lateral  pilot  model  was  converted  to  degrees  for 
proper  dimensionality.  Second,  the  output  force  commands, 

^LON’  ^LAT’  ^RUD*  were  given  in  units  of  pounds 

in  the  TAWDS  model,  were  assumed  to  be  multiplied  by  a gain 
of  one  degree/pound  before  being  fed  to  the  aircraft  flight 
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control  system. 

The  simulation  results  are  discussed  in  the  next  chapter. 
Tracking  error  time  histories  are  determined  against  targets 
in  straight  and  level  flight f and  JO,  60,  and  70  degree  bank 
turns.  In  addition,  the  effects  of  changing  the  pilot 
model  parameter  values  are  demonstrated. 


V.  Simulation  Results 


i 

I 
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Introduction 

In  Chapter  Twoi  the  basic  equations  of  the  Fire  Fly 
director  sunsi^ht  system  were  developed.  Chapter  Three 
presented  the  nonlinear  equations  of  motion  which  describe 
the  F-106b  aircraft.  In  Chapter  Four,  an  analytical  pilot 
model  for  air-to-air  aiunnery  was  developed  for  the  F-106. 

Thus,  the  necessary  elements  for  a closed  loop  dierital 
simulation  have  been  presented.  The  simulation  was  conducted 
usine:  the  EASY  Dynamic  Analysis  Programs.  The  closed  loop 
system,  including  the  equations  necessary  to  describe  the 
motion  of  the  target,  was  modeled  using  the  EASY  Model 
Generation  Program.  A complete  listing  of  the  closed  loop 
model,  and  a discussion  of  the  modeling  procedure,  can  be 
found  in  Appendix  A.  The  EASY  Analysis  Program  was  used  for 
the  time  simulation  of  manually  controlled  air-to-air  gunnery, 
the  results  of  which  are  presented  below. 

The  objectives  of  the  simulation  effort  were  to 
evaluate  the  terminal  air-to-air  gunnery  tracking  performance 
of  the  F-106  configured  with  the  Fire  Fly  manual  director 
gunsight,  and  to  adapt  the  TAWDS  multi-axis  pilot  model  for 
use  in  closed  loop  simulation  analysis.  In  addition,  the 
usefulness  of  the  EASY  programs  for  modeling  and  simulation 
of  dynamic,  nonlinear  closed  loop  systems  was  demonstrated. 

At  the  beginning  of  each  simulation  run,  the  attacking 
aircraft  was  trimmed  for  constant  altitude  flight  at  10,000 
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feet,  at  a speed  of  Mach  .8,  and  with  an  initial  bank  angle 
equal  to  that  of  the  target.  This  initial  trim  condition 
caused  the  stick  forces  at  the  beginning  of  each  simulation 
run  to  be  zero.  The  initial  position  of  the  target  with 
respect  to  the  attacker  was  set  so  that  the  elevation  and 
traverse  tracking  errors  were  near  zero.  The  initial  angle- 
off  of  the  target  with  respect  to  the  attacker  was  less  than 
two  degrees  for  each  simulation  run.  The  initial  target 
range  was  set  at  2000  feet  and  the  target  speed  remained 
constant  at  Mach  ,8.  The  target  maneuvers  consisted  of 
constant  altitude,  constant  bank  angle,  constant  rate  turns. 
Range  variations  for  each  10  second  simulation  were  less 
than  250  feet,  and  the  maximum  angle-off  encountered  was 
10  degrees. 

The  discussion  of  the  simulation  results  is  divided 
into  two  sections.  In  the  first  section,  the  results  of 
simulation  runs  against  targets  in  straight  and  level  flight, 
and  in  constant  rate  turns  of  30,  60,  and  70  degrees  of  bank 
are  discussed.  These  simulation  runs  were  conducted  using 
the  nominal  pilot  model  parameter  values  as  determined  in 
Chapter  Four.  In  the  second  section,  pilot  model  parameter 
values  are  varied  to  demonstrate  the  effects  of  these 
variations  on  tracking  performance. 

Simulation  results  are  presented  as  time  history  plots 
of  elevation  and  traverse  tracking  errors  and  the  corresponding 
longitudinal  and  lateral  pilot  stick  forces.  In  accordance 
with  the  conventions  specified  in  Chapter  Two,  positive 
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elevation  tracking  error  is  decreased  by  positive  pitch  rate; 
that  is,  by  positive  rotation  about  the  body-fixed  V axis,  as 
given  by  the  right  hand  rule.  Positive  traverse  tracking 
error  is  decreased  by  positive  yaw  rate;  that  is,  by  positive 
rotation  about  the  body-fixed  w axis,  as  given  by  the  right 
hand  rule.  Positive  longitudinal  stick  force  is  defined  as 
aft  stick,  the  direction  of  stick  movement  required  to  produce 
positive  pitch  rate.  Positive  lateral  stick  force  is  defined 
as  movement  of  the  stick  to  the  pilot's  right,  and  positive 
rudder  pedal  force  is  defined  as  right  rudder  pedal  deflection. 

Tracking  Response  with  Nominal  Pilot  Model  Parameter  Values 

In  this  section,  the  tracking  performance  of  the  F-106/ 
Fire  Fly  director  gunsight  is  discussed.  Target  maneuvers 
consisted  of  straight  and  level  flight,  and  constant  altitude, 
constant  rate  turns  at  30,  60,  and  70  degrees  of  bank.  The 
nominal  pilot  model  parameter  values,  as  determined  in 
Chapter  Four,  were  used  for  the  simulation  runs  presented  in 
this  section.  A listing  of  the  nominal  parameter  values  is 
presented  in  Table  II.  The  nominal  values  were  not  chosen 
to  obtain  the  best  tracking  performance,  but  were  selected 
to  be  within  the  range  of  values  recommended  in  the  TAWDS 
programs. 

The  discussion  in  Chapter  Two  indicated  that  for  the 
target  maneuvers  simulated  in  this  thesis,  the  time  required 
for  the  Kalmar,  filter  to  process  target  state  information 
has  a negligible  effect  on  tracking  performance. 
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Table  II.  Nominal  Pilot  Wodel  Parameters 


longitudinal  Kpg 

Parameters  * 

Kr 

''’e 

DZ. 

®LV 

Nominal  Values  1.0 

1.0 

1.0 

0.05 

5.0 

mr/sec 

0. 6,1,0 

Lateral  ^atl 

Parameters 

^RTR 

SuD 

''’t 

DZ* 

®LW 

Nominal  Values  O.36 

1.0 

0.0 

0.0 

0.05 

5.0 

mr/sec 

5.0 

deg 

This  fact  is  illustrated  by  comparing  the  tracking  error 
time  histories  for  a target  in  a 60  degree  bank  turn  with 
the  time  histories  for  the  same  target  maneuver,  but  with 
a simple  lag  transfer  function  with  a time  constant  of  one 
second  inserted  into  the  model  to  simulate  the  information 
processing  delay.  Figures  16  and  17  depict  the  tracking 
error  time  histories  with  no  delay,  and  Figures  18  and  19 
show  the  time  histories  with  the  one  second  delay.  Since  the 
error  time  histories  are  identical,  the  effects  of  the 
information  processing  delay  for  the  assumed  filter  model 
are  neglected  for  the  remainder  of  the  simulation  runs, 

A strict  judgement  as  to  whether  the  tracking  error  time 
“.•’'Tiee  presented  in  this  section  are  actually  similar  to 
*•  ••  - • jld  te  obtained  by  a human  pilot  flying  the 
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Fifijure  18.  Elevation  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains/1.0  Sec  Processing  Delay. 


same  aircraft/sight  configuration  against  the  sane  targets# 
cannot  be  made  until  actual  man-in- the-loop  simulation  data 
is  available  for  comparison.  However,  the  tracking  error 
time  histories  in  this  section  demonstrate  the  type  of 
tracking  response  predicted  by  McAir  (Ref  6i  66-69).  The 
elevation  error  response  contains  two  frequency  components. 

The  high  frequency  component  is  most  apparent  in  the  first 
four  seconds  of  each  simulation  run,  while  the  low  frequency, 
lightly  damped  component  is  clearly  present  throughout  the 
time  histories.  The  traverse  tracking  error  response 
consists  of  a single  frequency  component  which  has  a period 
of  five  to  nine  seconds.  The  mdst  significant  difference 
between  the  tracking  error  time  histories  obtained  in  this 
section  and  those  illustrated  in  Reference  6,  is  in  the  magni- 
tude of  the  tracking  errors,  particularly  the  traverse  track- 
ing error.  The  error  magnitudes  are  dependent,  at  least  to 
some  degree,  on  the  pilot  model  parameter  values,  as  will  be 
demonstrated  later  in  the  chapter. 

The  longitudinal  and  lateral  pilot  stick  force  time 
histories  depict  the  relative  stick  forces  applied  by  the 
pilot  model  where  the  reference  trim  stick  forces  in  both 
the  longitudinal  and  lateral  axes  are  zero. 

The  tracking  error  time  histories  and  corresponding 
stick  forces  for  a 60  degree  bank  target  were  shown  in 
Figures  16  and  1?.  Figures  20  and  21  present  plots  for  the 
same  target,  but  with  the  simulation  time  extended  to  20 
seconds  to  illustrate  the  limit  cycle  type  of  tracking  error 
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Fif^ure  20.  Elevation  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains/20  Sec  Simulation. 


del  Gains/20  Sec  SiiTiulation. 


response  over  the  longer  tracking  period.  Figures  22  and  23 
show  the  tracking  error  and  stick  force  responses  for  a 70 
degree  bank  target.  It  should  be  noted  that  the  scale  for 
plotting  the  longitudinal  stick  force  time  history  for  the 
70  degree  bank  target  in  Figure  22  is  doubled. 

Tracking  error  and  stick  force  time  histories  for  a 
30  degree  bank  target  are  shown  in  Figures  24  and  2$.  The 
traverse  tracking  error  response  depicted  in  Figure  25 
diverges  because  with  the  value  of  in  the  pilot  model 

set  to  1.0  as  in  Table  II,  the  lateral  clof  d loop  system  is 
unstable.  The  instability  can  be  corrected  by  reducing  the 
value  of  to  '^.5.  This  moves  the  closed  loop  zero  in 

the  lateral  pilot  model  closer  to  the  s-plane  origin  and 
thus  pulls  the  unstable  closed  loop  poles  back  into  the 
left  half  s-plane. 

For  a target  in  straight  and  level  flight,  the  F-106/ 
Fire  Fly  manual  director  gunsight  system  tracks  to  within 
one  milliradian  in  both  elevation  and  traverse  during  a 
10  second  simulation.  The  longitudinal  and  lateral  stick 
forces  are  negligible.  Thus,  plots  of  error  time  histories 
and  stick  forces  for  a straight  and  level  target  are  omitted, 

The  simulation  time  history  plots  presented  in  this 
section  indicate  that  the  tracking  responses  for  targets  in 
60  and  70  degrees  of  bank  are  stable  when  the  pilot  model 
parameter  values  are  set  to  the  nominal  values  shown  in 
Table  II,  Tracking  response  for  a target  in  30  degrees  of 
bank  is  stable  in  elevation  tracking  and  will  be  stable  in 
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Figure  22.  Elevation  Trackins;  Error  and  Stick  Force  vs.  Tine  for  70  De^  Bank  Target/ 
Nominal  Pilot  ^odol  Gains. 
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Figure  24.  Elevation  Tracking  Error  and  Stick  Force  vs.  Time  for  30  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains. 


Nominal  Pilot  Model  Gains. 


traverse  tracking  if  the  value  of  is  set  to  0.5.  These 

tracking  responses  must  be  compared  to  the  tracking  responses 
obtained  in  man-in- the-loop  simulations  with  the  same  aircraft/ 
sight  configuration  to  determine  their  true  validity,  but 
the  tracking  responses  obtained  here  do  exhibit  the  same 
structure  as  those  obtained  in  the  TAWDS  program.  In  the 
next  section,  the  nominal  pilot  model  parameter  values  are 
varied  in  simulation  runs  against  targets  in  60  degree  bank 
turns. 


Tracking  Error  Response  with  Varied  Pilot  Model  Parameters 

In  this  section  several  of  the  pilot  model  parameter 
values  are  changed  from  the  nominal  values  listed  in  Table  II, 
All  simulation  runs  are  conducted  against  targets  in  60 
degree  bank  turns.  The  longitudinal  and  lateral  pilot 
transmission  gains,  Kpg  and  are  varied.  Then,  the 

value  of  Kpyjj,  the  rudder  pedal  force  gain,  is  increased 
from  the  nominal  value  of  zero.  Next,  the  effects  of  adding 
differential  bank  angle  feedback  to  the  total  pilot  model  are 
determined.  Finally,  simulation  runs  are  conducted  with  a 
combination  of  different  pilot  model  parameter  values  to 
obtain  the  best  tracking  response.  No  claim  is  made  that  the 
responses  obtained  by  these  final  simulation  runs  are 
indicative  of  the  responses  which  can  be  expected  with  human 
pilots  flying  the  same  aircraft/sight  configuration.  The 
runs  were  made  strictly  to  show  the  effects  of  parameter 
value  changes  on  the  tracking  response. 
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ii. 

Figures  26  and  27  show  tracking  error  and  stick  force 
time  histories  for  a simulation  conducted  with  the  longitu- 
dinal and  lateral  pilot  transmission  gains,  Kp^  and 
set  to  the  minimum  values  for  which  the  tracking  response 
is  stable.  These  values,  as  determined  in  Chapter  Four,  are 
Kpg  *s  0.5  and  ~ 0*2.  The  plot  of  elevation  tracking 

error  in  Figure  26  shows  that  the  response  is  very  sluggish 
when  compared  to  the  corresponding  plot  for  the  nominal 
value  of  Kpg  = 1.0  in  Figure  l6.  It  is  noted,  also,  that  the 
lower  value  of  Kpg  does  not  excite  the  high  frequency 
component  of  the  elevation  tracking  error.  Comparing  the 
plot  of  traverse  tracking  error  in  Figure  27  to  the  corre- 
sponding plot  for  nominal  in  Figure  17  shows  only  minor 

differences  in  magnitude  and  period  of  the  response. 

Figures  28  and  29  show  tracking  error  and  stick  force 
time  histories  for  a simulation  run  conducted  with  Kpg  = 1.5 
and  = 0.54.  Comparing  the  plot  of  elevation  error  in 

Figure  28  with  the  plot  for  nominal  Kp^  in  Figure  16  shows 
a significant  reduction  in  tracking  error  magnitude  and 
a predominance  of  the  high  frequency  component  for  the 
higher  value  of  Kpg.  A comparison  of  the  traverse  tracking 
errors  in  Figures  29  and  17  shows  a modest  decrease  in  both 
magnitude  and  period  length  for  the  higher  value  of  K.,,. 

In  Figures  30  and  31»  the  tracking  error  and  stick  force 
time  histories  are  shown  for  a run  conducted  with  the  value 
of  Kpg  returned  to  its  nominal  value  and  further  increased 

to  a value  of  Kail-  ° .87.  In  this  case,  a comparison  of 
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Figure  26.  Elevation  Tracking  Error  and  Stick  Force  va.  Time  for  60  Deg  Bank  Target/ 
KpE  * 0.5  and 


Pifjure  29.  Traverse  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Kpg  = 1.5  and  - 0.54. 


traverse  tracking  error  responses  between  Figures  3I  and  1? 
shows  a marked  decrease  in  the  magnitude  of  the  error  and  a 
decrease  in  the  period  of  the  response  from  6 seconds  to 
approximately  4.7  seconds. 

In  Figures  32  and  33*  the  tracking  error  and  stick  force 
time  histories  are  plotted  for  a run  with  the  pilot  model 
parameters  returned  to  nominal  values  1 except  that  the  value 
of  Kj^ujj*  the  pilot’s  rudder  pedal  gain,  is  increased  from 
zero  to  Kj^ud  “ O.05.  Comparing  Figures  32  and  33  with 
Figures  I6  and  1?  shows  no  difference  in  the  tracking  error 
responses  and  a plot  of  rudder  pedal  force  commands 
in  Figure  34  shows  that  the  rudder  pedal  commands  follow  the 
pilot  lateral  force  commands  Figure  33*  t>vit  that 

the  rudder  commands  are  negligible.  This  corresponds  to  the 
claim  made  by  McAir  in  Reference  6 that  rudder  pedal  force 
commands,  when  used  at  all  by  the  pilot,  are  used  to  help 
coordinate  turns,  but  that  the  rudder  commands  are  generally 
negligible  in  the  terminal  phase  of  air-to-air  gunnery 
tracking.  Since  the  sideslip  angle  without  rudder  pedal 
force  commands  did  not  exceed  0,5  degrees  for  a typical 
simulation  run,  the  benefits  of  including  rudder  pedal 
commands  in  this  model  were  insignificant. 

Figures  35  and  36  depict  the  tracking  error  and  stick 
force  time  histories  for  a run  conducted  with  the  pilot 
model  parameters  set  at  the  nominal  values  except  that  the 
value  of  the  gain  on  the  differential  bank  angle  feedback 
was  set  to  ■ 0,1.  Comparing  the  elevation  tracking 


Elevation  Tracking;  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains/  Ko„p,  = O.O5. 
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Figure  33,  Traverse  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains/  = 0,05, 
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error  responses  of  Figures  35  ai-d  shows  a decrease  in 
the  magnitude  of  the  tracking  errors  for  the  response  in 
Figure  35*  Comparing  the  traverse  tracking  error  responses 
of  Figures  36  and  17  also  shows  a decrease  in  tracking  error 
magnitude  overall,  but  a significant  increase  in  the  initial 
time  required  to  achieve  zero  traverse  error.  The  changes 
in  the  tracking  error  responses  are  due  to  the  differential 
bank  angle  feedback  which  reduces  the  bank  angle  oscillations 
of  the  attacker  during  the  tracking  maneuver.  This  reduction 
in  bank  angle  oscillation  is  shown  by  comparing  Figure  37, 
where  the  attacker  bank  angle  is  plotted  versus  time  for 
the  simulation  run  with  differential  feedback,  to  Figure  38 
where  attacker  bank  angle  is  plotted  versus  time  for  the 
same  target,  but  without  bank  angle  feedback.  The  plots  in 
Figures  37  and  38  tend  to  support  the  theory  advanced  by 
McAir  in  the  TAWDS  program  that  the  pilot  uses  differential 
bank  angle  to  prevent  excessive  wing  rock.  However,  if  the 
pilot  model  uses  differential  bank  angle  feedback  when  the 
traverse  tracking  error  is  large,  the  result  is  an  increase 
in  the  time  required  for  the  pilot  model  to  drive  the  error 
to  zero.  Therefore,  if  the  pilot  uses  differential  bank 
angle  feedback,  he  probably  uses  it  only  when  the  traverse 
tracking  error  is  small. 

Next,  a simulation  run  was  conducted  which  combined 
relatively  high  values  of  Kp^  and  K^j^.  Kp^  was  set  to  1,5 

i 

and  K^jj_  to  1.0  to  show  the  effect  on  tracking  performance  j 

of  relatively  high  gains.  The  results  of  this  simulation 
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Fii^re  36.  Traverse  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Nominal  Pilot  Model  Gains  and  Differential  Bank  Angle  Feedback, 
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Figure  37.  Attacker  Bank  Angle  vs.  Time  for  60  Deg  Bank  Target  with  Differentia 
Bank  Angle  Feedback. 
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Figure  38,  Attacker  Bank  Angle  vs.  Time  for  60  Deg  Bank  Target  with  no  Differeni 
Bank  Angle  Feedback, 


run  are  shovm  in  Figures  39  and  40,  The  elevation  tracking 
error  response  depicted  in  Figure  39  shows  excellent  tracking 
for  the  entire  10  second  simulation.  The  traverse  tracking 
error  response  shown  in  Figure  40  shows  significant  improve- 
ment in  tracking  capability.  The  characteristics  of  the 
tracking  response  in  this  simulation  run  are  similar  to  those 
obtained  in  the  TAWDS  programs  except  that  the  traverse 
error  response  is  more  highly  damped  for  the  case  simulated 
herein. 

One  final  simulation  run  was  conducted  using  the  gain 

values  of  Kpg  = 1,5  and  = 1.0,  but  with  the  lateral 

pilot  model  deadzones  on  traverse  error  rate  (DZ«  ) and 

®LW 

differential  bank  angle  (DZ^^)  removed  from  the  model.  This 
technique  was  recommended  by  Wr.  Richard  Quinlivan  of 
General  Electric  Company,  Mr.  Quinlivan  is  the  General 
Electric  program  manager  for  the  Fire  Fly  program,  and  he 
contends  that  the  human  pilot  responses  do  not  exhibit  an 
appreciable  deadzone  on  either  traverse  tracking  error  rate 
or  on  differential  bank  angle.  The  results  of  this  simulation 
run  are  shown  in  Figures  41  and  42,  The  elevation  tracking 
error  response  shown  in  Figure  41  is  essentially  the  same  as 
the  response  obtained  in  the  previous  r\in.  The  traverse 
tracking  error  response  depicted  in  Figure  42  shows  a slight 
decrease  in  the  magnitude  of  the  tracking  errors,  but  the 
response  is  very  similar  to  that  obtained  in  Figure  40 
where  the  deadzones  were  not  removed.  Since  the  value  of 


in  this  simulation  run  was  zero,  removing  the 


deadzone 
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Fi(?ure  40.  Traverse  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  Bank  Target/ 
Kpg  = 1.5  and 


Figure  41.  Elevation  Tracking  Error  and  Stick  Force  vs.  Time  for  60  Deg  bank  Target/ 
^P£  “ ^*5  ^AIL  ~ ^*0/  Lateral  Pilot  Model  Deadzones  Removed, 
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on  differential  bank  angle  had  no  effect  on  the  traverse 
tracking  error  response.  However*  when  the  value  of  was 
increased  to  0.1,  the  tracking  response  observed  for  the 
traverse  error  was  essentially  the  same  as  the  response 
shown  in  Figure  36  where  the  deadzone  on  differential  bank 
angle  was  not  removed.  For  the  one  case  simulated,  the 
effect  on  the  traverse  tracking  error  response  of  removing  the 
lateral  pilot  model  deadzones  was  not  significant. 

Summary  of  Tracking  Performance 

This  chapter  presented  the  results  of  the  closed  loop 
simulation  of  the  F-106/Fire  Fly  manual  director  gunsight 
system  to  determine  the  tracking  performance  in  the  terminal 
phase  of  air-to-air  gunnery.  A multi-axis  pilot  model  was 
adapted  from  the  TAWDS  pilot  model  and  a nominal  set  of 
pilot  model  parameter  values  was  determined.  The  first 
section  of  this  chapter  discussed  the  results  of  simulation 
runs  against  targets  in  straight  and  level  flight,  and  in 
30,  60,  and  70  degree  bank  turns.  The  tracking  error  re- 
sponses were  found  to  have  characteristics  similar  to  those 
obtained  from  F-4  man-in- the-loop  simulation  studies,  except 
that  the  magnitudes  of  the  tracking  errors  observed  in  this 
simulation  study  were  larger. 

In  the  second  section  of  the  chapter,  simulation  results 
were  presented  from  runs  conducted  against  targets  in  60 
degree  bank  turns.  In  these  runs,  several  pilot  model 
parameter  values  were  varied  to  demonstrate  the  effects  on 
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the  tracking  error  time  histories.  The  values  of  the  longi 


tudinal  and  lateral  pilot  transmission  gains,  Kpg  and 
were  varied  from  the  minimum  values  necessary  for  stable 
tracking,  0.5  and  0,2  respectively,  to  values  of  Kpg  = I.5 
and  A significant  decrease  was  noted  in  the 

magnitudes  of  the  tracking  errors  as  the  gains  were  increased. 
The  addition  of  rudder  pedal  force  commands  in  the  lateral 
pilot  model  had  a negligible  effect  on  traverse  tracking 


error  response. 

The  addition  of  differential  bank  angle  feedback  to  the 
lateral  pilot  model  decreased  the  overall  magnitude  of  the 
traverse  tracking  errors  observed,  but  also  resulted  in  a 
marked  increase  in  the  initial  time  required  to  reach  zero 
traverse  error, 

A simulation  run  conducted  with  the  lateral  pilot  model 
deadzones  removed  showed  no  significant  difference  in  the 
characteristics  of  the  traverse  tracking  error  response. 

The  tracking  error  response  characteristics  were  shown 
to  be  highly  dependent  on  the  pilot  model  parameter  values. 
To  make  the  multi-axis  pilot  model  used  in  this  simulation 
more  useful,  the  parameter  values  should  be  adjusted  for 
good  correspondence  between  the  tracking  error  responses 
obtained  in  this  simulation  effort  and  the  tracking  error 
responses  obtained  from  F-106/Fire  Fly  man- in- the- loop 
simulation  studies. 

Table  III  presents  a set  of  recommended  pilot  model 
parameter  values  based  on  the  simulation  results  presented 


9? 


in  this  chapter 


Table  III.  Recommended  Pilot  Model  Parameter  Values 

for  F-106/Fire  Fly  Manual  Director  Gunsight 


Longitudinal  Kpp.  K- 

Parameters  ^ 

Kr 

Cry* 

®LV  ^ ^ 

Parameter  1.0-1. 5 1.0 

Values 

1.0 

0.05 

5.0  0.6,  1.0 

mr/sec 

Lateral  K.,, 

Parameters 

^RUD 

Tm  DZ*  DZ 

'lw 

Parameter  O.36-I.O  0. 5-1.0 

Values 

0.0 

0.0 

0.05  5.0  5.0 

mr/sec  deg 
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VI,  Conclusions  and  Recommendations  for  Further  Study 

Conclusions 

The  primary  objective  of  this  thesis  effort  was  to 
evaluate  the  tracking  performance,  during  the  terminal  phase 
of  air-to-air  gunnery,  of  the  Fire  Fly  manual  director 
gunsight  system.  Further  objectives  were  to  adapt  the  TAWDS 
multi-axis  pilot  model  for  use  in  the  F-106/Fire  Fly  simu- 
lation, and  to  determine  the  usefulness  of  the  EASY  computer 
programs  for  modeling  and  analysis  of  six  degree  of  freedom 
nonlinear  weapon  systems. 

To  accomplish  the  objectives  outlined  above,  a digital 
computer  model  of  the  aircraft/gunsight/ target/pilot  closed 
loop  system  was  developed  using  the  EASY  Model  Generation 
Program,  A set  of  nominal  pilot  model  parameter  values  was 
selected  from  within  the  range  of  values  recommended  for  the 
TAWDS  pilot  model.  The  tracking  performance  of  the  closed 
loop  system  was  simulated  using  the  EASY  Analysis  Program, 

Target  maneuvers  consisted  of  constant  altitude,  constant 
speed,  constant  rate  turns  at  bank  angles  of  30,  60,  and  70 
degrees.  After  establishing  that  the  system  model  could 
tracr.  the  targets  being  simulated,  the  pilot  model  parameter 
values  were  varied  to  note  the  effects  on  the  tracking 
error  responses. 

The  simulation  results  presented  in  Chapter  Five  show 
that  the  tracking  performance  of  the  Fire  Fly  manual  director 
gunsight  system  is  stable  for  the  target  maneuvers  simulated. 
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The  characteristics  of  the  tracking  error  responses  are 
significantly  dependent  on  the  pilot  model  parameter  values 
selected.  The  structure  of  the  tracking  error  responses  is 
similar  to  the  responses  obtained  for  the  F-4  aircraft  in 
the  TAWDS  simulation  studies. 

The  multi-axis  analytical  pilot  model  used  in  the 
simulation  is  readily  modeled  using  standard  components  of 
the  EASY  Model  Generation  Program,  If  the  pilot  model 
parameter  values  are  selected  to  obtain  a good  correspondence 
between  the  tracking  error  responses  presented  in  this 
simulation  study  and  the  responses  obtained  from  F-106/Fire 
Fly  man-in- the-loop  simulations,  the  pilot  model  should  prove 
to  be  a valuable  tool  for  future  simulation  studies. 

The  EASY  Model  Generation  and  Analysis  Programs  provide 
powerful  tools  for  modeling  and  analysis  of  complex  weapon 
systems.  The  EASY  Model  Generation  Program,  with  its 
capability  for  intermixing  EASY  standard  components  with  user 
Fortran  statements,  gives  the  programmer  the  necessax-y 
flexibility  to  accurately  model  a six  degree  of  freedom, 
nonlinear  aircraft  system  over  a wide  range  of  flight 
conditions.  The  EASY  Analysis  Program  provides  a wide 
variety  of  tools  for  the  design  and  evaluation  of  closed 
loop  systems.  The  ease  with  which  parameters  can  be  changed 
in  the  model  makes  the  program  particularly  useful  to  the 
control  systems  design  engineer. 
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Recommendations  for  Further  Study 


The  tracking  performance  of  the  Fire  Fly  manual  director 
gunsight  system  should  be  evaluated  against  dynarically 
maneuvering  targets.  This  task  can  be  accomplishea  in  two 
ways.  A target  subroutine  could  be  added  as  a standard 
component  of  the  EASY  Model  Generation  Program  and  then 
programmed  to  perform  the  type  of  maneuvers  desired  for 
simulation.  An  alternative  method  is  to  develop  a TAWDS 
model  of  the  closed  loop  system  and  run  simulations  against 
the  target  model  which  now  exists  in  the  TAWDS  program. 

For  simulations  against  dynamically  maneuvering  target 
models,  it  would  be  useful  to  explicitly  model  the  Kalman 
filter  and  the  tracking  system  of  the  Fire  Fly  system. 

Simulation  of  the  automatic  flight  control  capabilities 
of  the  Fire  Fly  system  could  be  performed  using  a model  of 
the  F-IO6,  but  an  appropriate  control  augmentation  system 
would  have  to  be  added.  Since  the  Fire  Fly  system  probably 
will  not  be  implemented  operationally  on  the  F-IO6,  it  would 
be  more  useful  to  model  an  aircraft  such  as  the  F-16,  which 
has  a control  augmentation  system  which  could  be  utilized 
for  analyzing  the  automatic  control  features  of  the  Fire  Fly 
system. 
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Appendix  A.  Computer  Modeling  and  Simulation  Techniques 

In  the  first  section  of  this  appendix,  a listing  of  the 

I closed  loop  simulation  model  is  presented.  Next,  the  func- 

tional flow  of  the  model  is  discussed.  Included  are  block 
diagrams  of  the  pilot  model/control  system  connections.  The 
I third  section  discusses  specific  techniques  used  in  the 

[ simulation  analysis. 
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Co*nputer  Wpdel 


cohhanos 


MOnet  OESCoioTICMsNONUIMEaR  F106  MODEL  «ITH  «»ILOT  MODEL ♦?!«£  CONTROL 
*00  T*«fLFS=OELCN9,-?3,';Y'’»,-E3tf^La’’,-23tCLNT,-23,CLOT,-25, 
CN*T,-2T,CMnT,-23,CNSA,-23,CNST,-23,C*OAa,-23,CNO*OE,-23, 
CY»T,-23,CrqT,-23 

*00  T«8LES3C70f6S.CLa,S4,CL<)A,68,CL0A2t*8|CL0B,Sa,CL0A,6ll, 
CLDff»68tCL'>t68,CLPA2t68tCLR>S8tCL<>A,<>8>CM0Et48,CM0,68fCN0,68, 
CHn»68fCNRA2f68tCNBR>63  » CNDA  »68 , CNnR,<}S  ,CNP  1 68 >CNR,68f  CYB,68 , 
CY9*2fS8tCYRB,68,CY0A,i 8,CY0R,6e,rYR,6S,CYR(68>CZA,6S,C20E>68, 
TN»«68|CLPA,6atC0RAGl,t52tC0RAG2>lR2,C3RAG3tl$2 
*00  P*RaME‘’ERS*  0LTANK,XCC,YREF,YC5,?REF,2CG,RPO 
*00  V*RIABLES=YAH,Cll,Cl2,C13,C21,C.?2,C23,C31,C32,C33,RAX,RAY,Ra7 
*00  VARIA3LFSaC3i  A,B,C,rO«-,PITCtROLL,aBV,aBH,RRY,Hg,aLPH,P,RH,n.VtELH 
*00  V*R1ARL^S=WRV,VRU,VRW,APU,ARV,APH,RTU,RTV,RTH,RA,RAD0T 
*00  V*RIAPLE5?=*TX,ATY,aTZ,vrx,VTY,VT-»,AAX,AAY,A*2,VAX,VAY,V*Z 
*no  V*RI*»L'SaARX,APY,*R2,VRX,VRY,7PZ,aAU,AAV,Aavi 
*00  7*R1A9LEE«*X,*Y,AZ 
L0C*TX0N>56  *V  tNPUT$«S0 

LOC*TION*76  FV  • INPUTS«$0< AL Y=X) , AV ( HAC=<» 

L0CATt0N*52  MA  INP0TS«FV (XsCl » , LGTT , 

L0C*TT0N»2  EN  INPUTS^MA (X*THR) 

FORTRAN  statements 
CCCC  •' 

CC  COMPOTE  ACCELEPATIONS  TN  G 

CCCC 


6a32.2 
*X>EII  AV/C 
*V>FV  *V/6 
*Z«EM  *y/G 

CCCC 

CC  COMPOTE  LONGITUDINAL  FOoCP  TERNS  FZ  EM  AMO  FX  EN 
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•>  CCCtf 

‘>  cccc 


INTERPOtftTE  TO  'TNO  LIFT  COEFFICIENTS 


> cc 
•>  cccc 
■» 

■> 

•> 

•» 

•> 

•>  cccc 
•>  cc 
•>  cccc 

•> 

•> 

•> 

•>  cccc 


CLOTO*  -TOLUl (MACAV,CL0T(4),CL0T(14I  , 1,-10) 

CL«TO«  -TOLUl  (MACAW, CLAT(4) ,CLAT(14) ,1,-10) 

CL0»  -TBLU2(MACA/,ALTSn,C70(9),C70t4) ,CZO ( 19) , 1 , 1 , -10 , -5 , 10,5) 
CLA«  -TeLU2(MACAW,auTSn,CZa(9),C?A(4)  ,CZa(l9),l,l,-10,-5,10,5) 
.5L0E»-T8LU2(MaCAW,aLTS0,CZ1E(9)  ,:ZOE  (4)  ,CZOE (19)  , 1 , 1 , -10 , -5, 10, 5) 

COMPOTE  LIFT,OPar,  COEFFICIENTS  ANO  SOOT  AXIS  FORCES 

ELE¥»  X2  LA  E’RPO 
ALPHA*  AL  AV*PPO 

CL1»  CLO  *•  CLOTO**  (CLA^CLATO) •ALPHA  • CLOE’FLEV 


•>  CC  LINEAR  INTERP.  0^  DRAG  COEFFICIENTS  AT  TNO  MACHS 

•>  CCCC 

•*  COOl*  T9LU2(X2  la  E,CL1,C0RAG2(28),C0RAG2(4)  ,C0P.AG2(33),1,1, 

•>  1 -9,-24,3,24) 

•*  IPCNACAV.LT.  .90)  GO  TO  100 

> :002a  TBLU2(X2  la  E,rLl,C3RAr;3(28),CDPaG3(4),C0RAG3(33),l,l, 

■>  t -5,-24,5,24) 

•*  C00»  (MACAW-. 90)»2a. *0002  ♦ ( . 95-MACA W) *20 . •COOl 

•>  in)  TO  200 

•>  100  CD02S  T9LU2(X2  LA  E, CL1,CQRAG1 (28 ) ,C0RAG1 (4) ,C0RAG1 (33) , 1, 1, 

•>  t -5,-24,5,24) 

> 500»  (.90-HACAV) •10.*CC02  ♦ (MACaW-. 80)*10 .’COOl 

•>  200  CALP*  COS(ALPHA) 

'>  SALP«  SIN(ALPHA) 

* PZSAW*  FZ  EN 

•>  FXSAW». FX  EN 

•*  *7  EHp-QS  AW*(CL1*CALP  ♦ CDO'SALi*)  ♦ FZ  EM 
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FX  Ell»  OS  AV*<CLl*SAUP  - C0tl»3ACP)  * FX  EN 


ecce 

CC  LOOK  UP  LONltTUOlNAL  OFPI'/ATIXES  FO?  P=Tr  AND  COMPUTE  TV  EM 

cccc 

MO  LO»  TBLU2«MACAV,ALTS0,Cf10(<»)  ,;*1014>  fCMOJlO)  ,1,1,-10,-5,10,5) 

MQ  L0a2.*T5LU2(MACAV,  ALTSO.rHIf  J)  ,010(4)  , CMOdO)  , t, t,- 10  ,-5,10  ,5) 
MAOLO«MQ  LO 

H0FL0»TBLU2(MACAV,lLTSn,CM0E(‘))  ,:H0E(4)  ,CMr)E  (19)  , 1 , 1, -10 , -5,10 ,5) 
CMOTO»  T9LU1  (MACA  V.CMOT (4 ) , ■;H')T(  14)  , 1 , -10) 

CHATO«  TBLU1(MACAV,CHAT(4)  ,'^HT  ( i4)  , 1 ,-l0 ) 

TV  EM*  OS  AV*C  l.O*(CMOTO  ♦ riATO* ALPHA)  ♦ TV  EN 
L0CATI0N*4  LO  IMPUTS=AV,LA  ^(X=ELE),EN 

FOMTBAM  STATEMENTS 
CCCC 

CC  A90  AEOO-FO?CE  TO’OUES  TO  H*TY2LO 

CCCC 


4 


TMP0*T9LU2(MACAV,ALTS0,TNP(9) ,TNP(4) ,TNP(19) ,1, 1 , -10, -5, 10 ,5) 
TV2LO*  TV2LO  ♦ (EREF-ZCG) •(^XELO-FXSAV)  - (TNPO*C  LO-XCG>* 

1 (FZ2LO-F2SAV) 


CCCC 

CC  COMPUTE  LATEPAL  OIoECTIOMAL  NONLINEAR  TERMS  FY1LD,TX1L0, TZILO 

CCCC 

CCCC  • 

CC  IMTEPPOLATE  TO  find  FY  DERIVATIVES 

CCCC 

5YBTO*  TBLU1(MACAV,CY3T(4) ,CV9T( 14),  1,-10) 

CYATO*  TBLUl(MACAV,CYAT(4),ryA-(14),l,-lO) 

CVPA20*  T9LU2(MACAV,ALTSO,CT9A?0),CYBA2(4),CY9A2(19)  ,1,1, 

1 -10,-3,10,5) 

5Y9AO*  TqLUl(MACAV,CV3A(4) ,CY9A ( 14)  , 1 , -10) 

CYBBOa  TBLU2(NACAV,ALTSC,CYB<»(P)  ,CY9B(4)  ,SVBq(19)  ,1,1, 

1 -10,-5,10,5) 

CC  COMPOTE  nonlinear  Y FORCES 
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Tir*1s  RLTANK’P  '50»RP0/(VT  AV  ♦ VT  AVI 
«TA  « 9E  AV*RPO 

FYILO*  QS  AV’CCYITO’RETA  ♦ CYATO’TERM  ♦ q€TA*(ALPHA* 

1 <CYBA20‘ALPHA  ♦ CYRAO)  » CYBBO* ABS ("ETA ») > 

ccce 

CC  IwrCRPOLATE  TO  PIMO  L»TX  DERIVATIVES 

CCCC 

CUPTO«  T9LUl(MACAV,r:LBT(4)  ,CL1T(14J 

CLBA?0»  T3tU2(lA0AVf  Al.TSn,CU9A2(9).CLB42C4)fCLBA2<19)  ,1,1, 

1 >10,-$, 10,5) 

:t9AO«T9LU2(MACAV,ALTSO,CLnA(9)  ,:u9A  (4)  ,01.  BA  (19)  , 1 , 1 , -10  ,-5,10  ,5) 
CLBBOa-'BLU2(MACAV,4LTSO,CL9B(<))  ,0199  (4)  ,01.93(19)  , 1 ,1,  - 10  ,-5, 10, 5) 
CtPA?0»  T8UU2(MAOAV,AtTSO,CLP12(9),OLPA2(4) , 0LPA2 ( 19) , 1, 1, 

1 -10,-5,10,5) 

CLPAOxTBLU2(MACA'/,ALTSO,CLPA  (9)  , OLPA  ( 4) , OLPA (19)  , 1 , 1 , -10 , -5, 10, 5) 
5LPAO*YBt.U2(MACAV,ALTSn,CLPA (9) ,SLPA (4) , CUP A (19) , 1 , 1 , -10 , -5,10 ,5) 

CCCC 

CC  COiPUTE  NOMLIVEAR  L=TX  TORaJES 

CCC5 

OLTAUK*  CC9T0*BETA  -(CL0T0*9LTAN<  ♦ CYATO»2 .33) ’TERM 
CNOLIHs  ((CLBA20»A9S (ALPHA)  ♦ rL9AO) •ALPHA  ♦ CL9BO*ABS (BETA) ) 

1 •BETA  ♦ ((CLPA20*ALPHA  ♦ CLPAO)*’  SO  ♦ CLRAO*R  SO) 

? •ALPHA’S  LO*RPO/(VT  AV*VT  AV) 

TXILO  * OS  AV*3  L0*(CL'^ANi<  ♦ CNOLIM) 

CCCC  '■,  * 

CC  INTERPOLATE  T3  pTND  N»TZ  DERIVATIVES 

CCCC  , 

CN9TO»  T9LU1(MACAV,CN9T(4)  ,0N9’‘ ( 14)  , i ,-10) 

CNBA20«  TBLU2(MA0AV,ALTS0,CNPA?(9),CN9A2(4),CM9A2(19),1,1, 

1 -10,-5,10,5) 

CNRAO*  T9LU1(KACAV,CN9A(4) ,CN9A(14) ,1,-lC) 

CNB9(’«YBLU2(MACAV,ALTS9, 0599(9) ,;n1B ( 4) ,0599 (19) , 1 , 1, -10 ,-5, 1 0 ,5) 
OFLCMO*  T9LU1(HA0AV,9ELCN9(4) ,OEL0N9(14) ,1,-10) 
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CNOAAOs  TBLUl  CNnAA((t)  .rMOAAdA) 

COAO^Oa  TBLUl (HACAV , CNOAOB (4) ,r NBAOE ( 14)  . 1 1 -10) 


> CCCC 

* CC.  C01«**JTE  HONLINPAI?  MaTZ  TOBQUES 

> CCCC  , 

> CHTAtUta  CNBTO’OETA  - 0.6*TEBH 

> ATL«  X2  LA  A»RPO 

* CNONLIs  ALPHA’COBAZO^ALPHa  ♦ CN3AO)  ♦ BET  A*  (CN89O*A0S  CBET  A) 

>1  ♦ OELCNO)  ♦ (CNnaAO*AL!»HA  ♦ COAOEO^ELEV) ‘AIL 

* TZIL3»  os  4V»B  LO’tCMTftNX  CNONLI) 

> CCCC 

> CC  LOOK  UP  LATE RAL-OI^ECT TONAL  LINEAR  COfeFEICIFMTS  FOR  LO 

> CCCC 

> VB  L0aT8LU2(HACAV,ALTSn,CYO(R) ,CTBt4) ,CYB( 19) , 1 , 1 , -10 , -5, 10,5) 

> Y»  L0aT9LU2»HACAV,aLTSn,CYP(9),CYP<4) ,CYO(19),1,1,-10,-5,10,5) 

> YR  LOaT8LU2(NACA7,ALTSn,CYP(?),CYP<4) ,CYR(19),1,1,-10,-5,10,5) 

> Y*)aLnaT8LU2(M'aCAV,ALTSO,CYnA(9),:YnA{4),CY0a(l9)  ,1 , 1, -10 , -5, 1 0 , 5) 

> YnRL''aTBLU2(MACAY,ALTS'»,CY')R(9),CY0R(4),CY0R(19)  , 1 ,1 , -10  , -5,  10  , S) 

* > LP  LOaTBLUZtMACAYjALTSn.CLotRl.CLP***)  ,CLP<19),1,1,-10,-5,10,5) 

I 

j > LR  LnaTBLU2(HACAV,ALTS0,CLR('>),CLR('»)  ,CL?(19),1,1,-10,-5,10,5) 

\ > lb  LO»TeLU2(MACAV,ALTSn,rLB(P),CLB{4) ,CL9(19),1,1,-10,-5,10,5) 

I 

I > LOAL'’aTBLU2(PACA'/,AL''SO,CLOA(9)  , CLOA  ( 4)  ,CLOA  (19)  , 1 , 1 , -10  , -5,10 , 5) 

» . L0RL6aTBLU2(»1ACA7,ALTS0,:LPR('^),:L0R(4),CL0R(19)  , 1 , 1, - 10 ,-5,1 0 ,5) 

* NT  LOaT8LU2{MACA\(,ALTSO,CNPl9), 099(4)  ,CN3  ( 19 ) , 1 , 1 , -10  , -5 , 1 0, 5) 

* >»OALOaTeLU2(MACA7,ALTSO,CNOA (9) ,CNOA (4) , CUBA (19) , 1 ,1, -10 , -5 , 10 ,5) 

» NnRLnaTBLU2(MACAV,ALTSO,CNOR (9) ,:nOR(4) ,CN0R(19) , 1 , 1, -10 ,-5 , 10 ,5) 

> RP  L0aTBLU2(MACAY,ALTS0,CNP(9),CNP('»)  ,CNP  ( 19)  , 1 , 1 , -10  , -5 , 10,5) 

> MR  LnaTBLU2(MACAY,ALTSO,CNR(9),ONP(4) ,rNR(19),l,l,-10,-5,10,5) 

> LOCATIONa2B  LO  INPUTSaAV,LO,L A A(XaAIL)yLA  R(XaRUO) 

> FORTRAN  statements 

> CCCC 

> CC  AOO  AERO-FORCE  TOROUES  TO  L*TX2LO  ANO  NaTZZLO 

> CCCC  • 


7 
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TX2tO»  TX2tO  ♦ (YREF-VCr.J»(rz?|.0-FZSAW)  - ( ZREF-ZCC) •FrZLO 
TZ2L0»  TZZtO  ♦(TNPO*C  LO-XCG)* ^TZLO  -(VRPP-YCG) • <FX2L0-FXSAV» 
L0aATI0»|s9  SO  INPUTSsLOtLO 

LOCATIOM»AO’*AFAX 

l,OCATIO>la40«‘,UAAX,INPUT3*AFAX  * 

LOCATIO>l>AO<i>AFAY 

UOCA TIOSaAO ^,L AAY, inputs »AFAY 

l0CATI0MaA2J,AFAZ 

».0CATI0M*42-,LAAZ,IMPUTSaAFAZ 

LOCATIONaAZOf AFtfX 

L0CATI0N344C,AFVY 

L0CATt0N*4AC,AFVZ 

FORTRAN  STATEMENTS 

CCCC 

C COMPUTF  TARGET  VELOCITIES  ANO  ACCEL  IN  XfY.Z 

CCCC 

ATXaX?  LAAX 
ATV»X2  LAAY 
ATZ»X2  LAAZ 
WTX*X1  AFVX 
VTYaXl  AFVY 
VTraXl  AFVZ 

fees . 

C transformation  MATRICES 

CCCC  » 

YAMaYAHSO*R°0 

PITCboitsO*RPO 

ROLLaR(X.SO*RPO 

CllaCnSCYAMjaCOSi^ITC) 

C12aSIN(YAH»*C0S(»ITC) 

CH«-RIN(PITC) 

CZlaCOSCYAMI^SINCaiTO'SINtPOLU  -SIN(  YAM) ’COStROLL) 
C22»STM«YAM)»S1N«“ITC) ’SIN (POLL) ♦COS (YAM) 'COS (ROLL) 
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' C?3»COS{PITC»*SIN<ROLL) 

C31«C0'! C YAM) ♦SIN(=»ITC) ’COS  (ROLL)  ♦SIM{ YAM)  •SIM  (ROLL ) 
C32»SIN ( YAM) •SIN(3ITC) ‘COS (ROLL) -;0S( YAM) ’SIM (ROLL) 
C33«COS(PITC)*COS(ROLL) 

CCC8C 

C COHPOTF  ATTACKER  ACCELERATIONS  IN  U,ViM',AXES 
CCCCC 

AAUmU"  LO-EU  AV-32.a*SIN(PITC) 

AtVaVO  LO-EV  A V«3? .2»COS (PITC) 'SIN (ROLL) 

AAMsMO  LO-EH  AV>32.2*COS(PITC) ‘COS (ROLL) 

ccce 

C COMPUTE  ATTACKER  YELOCITIES  AMO  ACCEL  IN  X,Y,7 
CCCS 

AAXaCll»AAU+C21*AAV»C3l*AAM 

AAY»C12*AAU*C22*AAV»C32»AAM 

AA7»C13»AAU»C23*AAV»C33»AAM 


VAX«Ctl*U 

SO*C21*V 

SnfC31*M 

SO 

WAY»C12»U 

SO ♦CEE* V 

S0*C32*M 

so 

VAY«C1T*U 

SO*C23*V 

S0*C33*M 

so 

CCCC 

C COMPUTE  RELATIVE  ACCEL  AND  VELOCITIES  IN  X.Y.r 
CCCC 

•ARX«ATX-AAX 

ARY»A’'Y-AAV 

A«»A*’7-AAZ 

¥RX«VTX-VAX 

VRY«V^Y-VAY 

VR7»VTZ-VA7 

CCCC  • 

C CONPU'E  RELATIVE  ACCEL  AND  VELOCITIES  IN  ROOT  AXES 
CCCC 

ARU»Cll*ARX*Cl2*ARY»C13»ARr 

ARV«C21^ARX*C22*ARY»C23»AR7 
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kW^CZi.* 


VRU»C11*VRX»C12*V^Y*C13»V'?7 

VRV»C?1»VRX»C22*WRY»C23»VRY 

VRW«C31»VRX*C32*VRY*C33»VRY 

CCCC 

c co*tPuTe  range  rate 
CCC5 

»A')3T=S0RT(YTX**2* VTY**2+VT7**2) -SORT ( YAX**2*VAY* *2*VAZ** 2) 
XI  ITRXaVRX 
XI  ITRYaVRY 
XI  tTR2*VR2 
l,OCAT10N««*3ft,ITRX 
L0CATt0N»i.3^,MARX,IMPUTS*ITRX 
tO<;ATI0Na45S,ITRY 
L0CATI0N«45H,'IARY,INPUT5aITRY 
t0CATI0N«47G,ITRZ 
LOCATION»47A,MARZ,INPUTSsITR2 
PORT RAM  statements 
CCCC 

C COMPUTE  RANGE  COMPONENTS  IN  PIXEO  COORDINATES 
CCCC 

RAXaXZ  MARX 
r<aY«X2  MARY 
RAZ«X2  MARZ 

CCCC 

c compute  range  magnitude 

CCCC 

•AaSDPTimx^^Z^RAYa’Z^RA^aaZ) 

CCCC 

C TRAHSPORM  RANGE  COMPONENTS  TO  SOOT  AXES 

CCCC 

R^tl«Cll»RAX»C12*RAY»C13*RAZ 

RTV«C21*<^AX*C22*RAY»C23*RAZ 
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* RTW»C:»1*HAX»C32*R4Y*C33»RA7 
» CCCC 

> C .COMPUTE  BULLET  TIME  OF  FLIGHT 

> CCCC 

> C^«C.00625*SIGAV»$ORT(VT  4V*33G0.) 

> *«CP*(VT  AV  . A DOT- 13. 4) 

» B»CB*OA>RAOOT-3341 . 

» C»RA 

> T0F»(-9-SQRT(B**2.-4.*A*C) )/(2.*AI 

> CCC5 

> C COMPUT  RANGE  COMPONENTS  REOUIPED  FOR  HIT 

> CCCC 

> ABV»VO  LO-EV  AV 

> '»PV»-”OP*(VRV>0.5*TOF* ( ABV+APV)) 

> HB»l/(l»CP*TOF) 

> ABMsHD  LO-EH  AV 

> ALPH*AL  AV*RPO 

> RRH«-“OF*(H9*3350.*(-. 04362)  ♦VEH*0 .5*T0F* ( A9H>ARH) 

> 1 ♦T0F“H9*'CB*VT  AV*  SIN  ( ALPM) ) 

> CCCC 

» C COMOUTF  ERROR  SIGNALS  TO  ORIVF  “ILOT  MODEL 

> CCCC 

> ,ELV«(RRH-RTH)/RA 

> ELN»-(R?V-RTV) /RA 

* C2,  MAPPsELV 

> C2  MAOO«ELM 

* L0CATI0N*542,MAPo,INPUTS»ITLP (*=X) 

* L0CATI0N*934,MALLtINPUT5s'1APP 

> L0CATI0N*536,ITLPtINPUTSsMALL 

> LOCATIONaSS'fTFLPtlNPUTSalTLP 

> locatiom»fip,salp,inputs=mall 

> LOCATIOMaSTOtMALPtlNPUTSsITLP.SALPCXrCE) 

* IOC6TIOM»560,‘1CLP,INPUTS»MALP,TFLP 


I 
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> LO';»TtO‘l»5S3,TTAI,lMPUT3  = MCLP 

»>  U0C»TI0Na6l.l,MAn0,INPUT3  = ITL0(X=XJ 
•>  LOCA TI0N»6T'f  ■< ACL,  INPUTS sMA 00 
•>  LOCATrON=63':,rTLO,INPUTS=MAOL 
•>  L0CATI0Ma6L'=, SA03, INPUTS  =MA0L 
•>  L0CATI0M=63', HCOO, INPUTS =SA 00, ITLO 
■*  L0CATI0N=e02,»<CPA,INOUTS=S0<R0L=X) 
•>  LOCATION=60t,SALO,INPUTS=MCOA 
■*  LOCATION=639,MCLO,INPUTS=SALO,MCOO 
■>  COCA TI0N=65',»1A02, INPUTS  =1000 

> L0CATI0N=6SO,NA01,lNPUTS=MCt0 


tOCATION=102 

FUTQ 

lNPUTS  = Al/{SIG  = X) 

tOCATIONsli.2 

FUXQ  . 

INPUTSaAVf  (STG  = XI 

L0CATI0N=1?2 

LE  Q 

INPUTS=S0(0=X) ,FUTQ(X=PO) ,puKQ(X=GAI) 

L0CATI0N=12t 

SA  Q 

INPUtSaLE  0, 

LOCATION=lAfc 

LA  0 . 

INPUTS=SA  0, 

tOCaTlO»=16«.,MC 

E, INPUTS 

=LA  0,ITAI,LSET(X=C4> 

L0CATI0M=166 

SA  E 

INPUTSaMC  E, 

L0CATI0N=1&8 

LA  E 

INPUTS=SA  E, 

L0CATI0N=223 

FU  A 

INPUTSaAW(OC=X) 

L0CATI0Na2C 1,M A AI, inputs 

=LGAT(X=C2),‘1A01{XsX) 

L0CATI0N=2C3,MA 

A, INPUTS 

= FU  A (X=r;i)  ,NAAI 

tOCATION=20r 

LA  A 

. INPUTS=MA  A, 

L0CATI0N=233 

LE  R 

IMPUTS=SO(’=X) 

L0CATI(VI=?S3 

FU  P 

INPUTS=AV(0:=<) 

t0CATI0N=2S5 

LE  P 

INPUTS=SO(»=X) ,FU  PtX=GAI) 

L0CATI0N=235 

HC  R 

INPUTSaHA  A,LE  R,LE  P,LGRT(X=C4) 

U0CATI0N»237,NA  R, INPUTSslA 02 , NC  P(X=S2» 
I.0CATI0N»239,LA  i^tlNPUTSaiA  R 
L0CATI0N3317,MAAT, INPUTS =S0 (ROL^X) 

LOCA TIONa332,LCAT, INPUTS *iA AT, 
L0CATI0M=31A,'*ART,INPUTS=AW  (EV»X» 
to: ATI0N*3»V,tGRT, INPUTSaH ART, 
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> LOC»TrONs71’,MAET,INPUTS=SO  (ALTsX) 


* L0CATI0M=3’'^,LGET,INPUTS=*1AET, 

> LO5ATION=3?0,MATT,INPUTS=AV  (lACa*) 

* LOCATrON=3-0,LGTT,INPUTS=HATT, 

> EMO  OP  MOP'’L 


> P^XMT 
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Discussion  of  the  Cowputer  Model 


I 


The  EASY  computer  source  program  for  the  closed  loop 
F-106/Fire  Fly  gunsight  system  was  presented  above.  This 
source  program  contains  standard  components  of  the  EASY  Model 
Generation  Program  intermixed  with  Fortran  statements.  A 
listing  of  the  EASY  standard  component  names  used  in  this 
program  and  their  description  is  presented  in  Table  IV. 


Table  IV.  EASY 

Standard  Components  (from  Ref  12 i 22-23) 

Standard  Component 
Name 

Description 

AF 

Analytic  Function  Generator 

AV 

Aerodynamic  Variables  from  States 

EN 

Engine  Model 

FU,  FV 

Tabular  Function  Generators,  One 
Input,  Two  Inputs 

IT 

Integrator 

LA 

Lag  Transfer  Function  (Time  Constant 
Form) 

LD 

Lateral  Aerodynamic  Model 

LE 

Lead  Lag  Transfer  Function  (Pole- 
Zero  Form) 

LG 

Lag  Transfer  Function  (Pole  Form) 

LO 

Longitudinal  Aerodynamic  Model 

MA 

Multiply  and  Add  (2  inputs) 

MC 

Multiply  and  Add  (4  inputs) 

SA 

Saturation  Function 

SD 

Six  Degree  of  Freedom  Rigid  Body 
Dynamics 

TF 

Transfer  Function  (First  Order 
Numerator  - Second  Order  Denominator) 
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In  naming  EASY  standard  components  and  the  physical  quantities 
associated  with  each  component,  the  following  character 
assignments  are  specified  (Ref  12i  15) « 

7 Character 
Name 

Physical  Standard  Specific 

Quantities  Component  Component 

Name  Identifier 

1 

Specific  Component  Name 

Thus,  a lag  transfer  function  in  pole  form,  LG,  could  be 
identified  simply  as  LG.  If  more  than  one  lag  transfer 
function  is  needed  in  the  model,  separate  lag  transfer  functions 
would  be  identified  using  the  character  specific  component 
identifier.  Hence,  LGET  identifies  a lag  transfer  function 
used  for  elevator  trim,  and  LGRT  identifies  a lag  used  for 
rudder  trim.  The  three  character  physical  quantity  name  aids 
in  making  connections  between  components  and  in  specifying 
parameter  values.  For  example,  ALTSD  is  the  output  state 
which  identifies  the  aircraft  altitude  as  computed  in  the 
six  degree  of  freedom  rigid  body  dynamics  standard  component. 
GAILAAX  identifies  the  gain  parameter  for  the  lag  transfer 
function  in  time  constant  form,  LAAX, 

The  use  of  EASY  standard  components  allows  the  programmer 
to  model  most  of  the  basic  aircraft  equations  and  flight 
control  system  components  simply  by  specifying  input  and 
output  connections  for  the  components  and  by  specifying 


T » 
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r 

the  component  parameters.  The  use  of  Fortran  statements 
allows  flexibility  in  the  use  of  table  look-up  functions 
for  aero-coefficients  and  in  the  specification  of  nonlinear 
stability  derivatives. 

The  EASY  model  program  used  in  the  simulation  is  divided 
into  sections  as  indicated  by  the  numbers  along  the  right 
hand  side  of  the  program  listing.  In  Section  One,  the  basic 
connections  between  the  aerodynamic  variables  component  (AV) 
and  the  engine  model  component  (EN)  are  specified.  The 
component  FV  is  a two  dimensional  table  look-up  function 
which  is  used  to  determine  a nominal  thrust  level  as  a function 
of  altitude  (ALTSD),  and  Mach  number  (MACAV).  This  nominal 
thrust  value  is  multiplied  (MA)  by  an  incremental  thrust 
value  which  is  determined  by  the  thrust  trim  component  LGTT. 

The  resulting  required  thrust  level  (THREN)  is  fed  into  the 
engine  component  (Ref  12 i 117). 

Section  Two,  written  in  Fortran,  provides  interpolation 
functions  to  determine  the  lift  coefficients  from  tables  of 
data.  The  general  form  of  the  table  interpolation  function 
can  be  found  in  Reference  12.  Then  the  equation  to  determine 
the  total  lift  coefficient,  CLl,  is  specified. 

In  Section  Three  the  drag  coefficients  are  determined 
by  interpolation  and  the  total  drag(CDO)  is  determined.  Then 
the  total  longitudinal  force  equations  (FZ  EN  and  FX  EN)  are 
specified. 

In  Section  Four,  the  pitching  moment  coefficients  (C  , 

o 

^m  * ^m  * ^m  ^ determined  by  interpolation  of  tabular 

I 1 

I 

i 

L 
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values  and  the  tank  pitching  derivatives  (C  and  G ) are 

m 

added  to  the  engine  pitching  moment  to  yield  the  total 
external  moment  TYEN.  Next,  the  standard  component  for  the 
longitudinal  aerodynamic  model  (LO)  is  specified  which 
determines  the  linear  pitching  moment  equation.  The  aircraft 
neutral  point  (TNP)  is  found  by  interpolation  and  used  in  the 
equation  for  the  total  pitching  moment  TY2L0  (Ref  12 i 119). 

In  Section  Five,  the  nonlinear  side  force  derivatives 


(C  , C , C 
Yfl  y. 

Ht 


0-m  ^3  2 

T a 


, C , and  C 
^3„  ^3  l3( 


are  determined  by 


interpolation  and  used  to  determine  the  nonlinear  side  force 
FYILD.  In  Section  Six,  the  nonlinear  rolling  moment  derivatives 


(C 


2 * ^2 


, C,  , C.  , C.  , C.  , and  C.  ) are  found 
3^  3 |3l  P^2  ‘‘P.  ‘‘r. 


T ■'a”  "a  ^ "a'  "a  ”a 

by  interpolation  and  used  to  determine  the  nonlinear  rolling 
moment  TXILD. 

Section  Seven  determines  the  nonlinear  yawing  moment 


derivatives  (C  , C 
nQ  n 
Pi 


3 2 


. , C, 


n 


3, 


""3  131 


. AC  , C. 
^3 


, and  C 


n 


) 


T ^a-  ■'a 

by  interpolation  and  uses  the  values  to  determine  the  nonlinear 

yawing  moment  TZILD.  In  Section  Eight,  the  linear  lateral 

side  force  coefficients  (C„  t , C„  , and  C„  ),  the 

^3  ^P  ^R  ^6^^ 

linear  rolling  moment  coefficients  (C,  , C,  , C,  , C,  , and 

*P  R 

C,  ),  and  the  linear  yawing  moment  coefficients  (C  » C , 

Hq  n. 


Cn  • C , and  C ) are  determined  by  interpolation  and  then 

lie  llT^  n«s 
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the  lateral  directional  aerodynamic  model  standard  component 
(LD)  is  specified.  In  Section  Nine*  the  total  rolling  and 
yawing  moments  TX2LD  and  TZ2LD  are  determined  by  adding  the 

1 

I 

moments  due  to  the  eg  moment  arms.  Then*  to  complete  the 
basic  aircraft  model*  the  standard  component  SD  is  specified 
to  determine  the  equations  of  the  six  degree  of  freedom  rigid 
body  dynamics. 

Section  10  specifies  the  target  acceleration  and  velocity 
components  in  inertially- fixed  coordinates  XY2.  The  AF 
functions  specify  analytic  equations  for  the  acceleration 
and  velocity  of  the  target.  For  the  simulation*  the  target 
maneuvers  are  constant  speed*  constant  altitude*  constant 
rate  turns.  Thus*  the  target  acceleration  vector  is  given  by 

» -A  8in((tf^t)  X + A cos(u^t)  X (75) 

where  A is  the  magnitude  of  the  target  acceleration  and 
is  the  target  turn  rate.  The  fixed  coordinate  system  xXZ 
has  its  origin  at  the  center  of  mass  of  the  attacker  and  is 
directed  so  that  at  time  t > 0*  X is  out  the  nose  of  the 
attacker*  Y is  out  the  right  wing*  and  2 is  determined  by  the 
right  hand  rule. 

Similarly*  the  target  velocity  vector  is  given  by 

• V cos(tt)^t)  X + V 8in((i)^t)  Y (76) 

where  V is  the  magnitude  of  the  target  velocity. 
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The  components  AFAX,  AFAY,  and  AFAZ  are  the  X,  Y,  and  Z 
components  of  the  target  acceleration  vector.  The  components 
LAAX,  LAAY,  and  LAAZ  are  simple  lag  transfer  functions  which 
£ ulate  the  time  delay  in  the  processing  of  target  acceler- 
a>.  ^n  information  by  the  Kalman  filter  in  the  Fire  Fly  system. 
The  components  AFVX,  AFVY,  and  AFVZ  are  the  Xt  Yt  and  Z 
components  of  the  target  velocity  vector.  The  components  of 
target  acceleration  and  velocity  are  then  renamed  for 
convenience. 

In  Section  lit  the  components  of  the  transformation 
matrix  to  transform  between  fixed  coordinates  XYZ  and  aircraft 
body  coordinates  uw  are  specified  in  terms  of  the  Euler  angles, 
YAW,  PITC,  and  ROLL.  The  parameter  RPD  converts  the  Euler 
angles  from  degrees  to  radians.  In  Section  12,  the  attacker 
acceleration  components  are  computed  in  body-fixed  coordinates 
and  transformed  to  fixed  coordinates.  The  attacker  velocity 
components  in  body-fixed  coordinates  are  also  transformed  to 
fixed  coordinates.  Then  the  components  of  relative  acceler- 
ation (ARX,  ARY,  and  ARZ)  and  the  components  of  relative 
velocity  (VRX,  VRY,  and  VRZ)  are  computed  and  transformed 
back  to  body  axes. 

In  Section  13i  the  range  rate  (RADOT)  is  computed.  The 
equation  for  RADOT  is  strictly  valid  only  if  the  attacker  and 
target  velocity  vectors  are  collinear.  In  the  terminal  phase 
of  tracking  considered  here,  the  velocity  vectors  remain 
essentially  collinear,  and  thus  the  given  equation  is  used. 

The  relative  velocity  components  in  fixed  coordinates  are 
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integrated  and  added  to  the  initial  range  components  (specified 


in  MARX,  MARY,  and  MARZ)  give  the  present  range  in  fixed 
coordinates.  The  presen^  range  components  are  used  to  compute 
the  present  range  magnitude  and  are  then  transformed  to  yield 
the  present  range  components  in  body- fixed  coordinates  (RTU, 
RTV,  and  RTW). 

Section  14  computes  the  bullet  time-of- flight  (TOF)  and 
the  range  components  required  for  a hit  (RRV  and  RRW).  The 
required  range  components  are  compared  to  the  present  range 
components  to  yield  the  elevation  and  traverse  error  signals 
to  be  fed  to  the  pilot  model. 

Section  I5  specifies  the  multi-axis  pilot  model  in 
terms  of  EASY  standard  components.  Figures  43  and  44  are 
block  diagrams  which  show  the  connections  of  the  longitudinal 
and  lateral  pilot  models  with  the  longitudinal  and  lateral 
flight  control  systems.  In  parentheses  near  each  component 
are  the  EASY  acronyms  which  represent  each  component  in  the 
model.  The  parameter  values  shown  are  the  nominal  values 
which  were  used  in  the  simulation. 

The  standard  components  in  Section  1?  are  used  for 
trimming  the  aircraft  at  the  beginning  of  each  simulation 
run.  After  an  initial  trim  condition  is  reached,  the  output 
states  of  the  lag  components  for  aileron  trim  (LGAT),  rudder 
trim  (LGRT),  elevator  trim  (LGET),  and  thrust  trim  (LGTT) 
are  held  constant.  These  trim  components  are  included  in 
Figures  43  and  44  and  are  outlined  by  dashed  lines. 
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Discussion  of  Simulation  Techniques 

This  section  presents  some  of  the  specific  techniques 
used  in  the  simulation  program.  For  more  complete  informa- 
tion on  the  use  of  the  EASY  Analysis  Program  the  reader  is 
referred  to  References  12  and  I5.  Included  in  this  section 
is  an  outline  of  the  procedures  used  for  a typical  simulation 
run,  and  a discussion  of  some  of  the  sensitive  parameters  in 
the  model  and  how  the  parameter  values  can  affect  computer 
run  costs. 

A typical  simulation  run  begins  with  the  specification 
of  all  the  parameter  values  and  tables.  Next,  initial  condi- 
tions are  specified  to  define  an  operating  point.  Typically, 
these  initial  conditions  include  initial  altitude,  airspeed, 
and  bank  angle.  Next,  the  steady  state  command  is  used  to 
specify  values  of  the  states  at  the  desired  trim  point.  The 
operating  point  state  values  are  then  loaded  into  an  initial 
condition  vector  for  iiiitializing  a time  simulation.  ! 

The  computer  simulation  time  is  very  sensitive  to  the 
arrangement  of  the  model  components.  Derivative  operations 
in  the  forward  loop,  such  as  in  the  pilot  model,  cause 
computer  run  time  for  a ten  second  real  time  simulation  to 
be  near  400  octal  seconds.  By  eliminating  the  specific  deri- 
vative operations,  simulation  computer  run  time  is  decreased  j 

by  almost  a factor  of  ten.  In  the  model  used  for  the  simula- 
tions, the  pilot  model  requires  as  inputs  the  derivatives  of 
the  tracking  errors  in  elevation  and  traverse.  Rather  than 
trying  to  specifically  take  the  derivatives  of  the  errors. 
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simple  lag  feedback  circuits  were  modeled  and  the  error  rates 
were  picked  off  just  ahead  of  the  integrators.  This  concept 
is  illustrated  in  Figures  43  and  44.  Care  must  be  exercised 
in  the  selection  of  the  time  constants  used  in  the  feedback 
circuits.  These  time  constants  (Tg  and  T^)  represent  the 
pilot's  lag  in  response  to  the  observed  elevation  and  traverse 
tracking  errors.  Making  the  values  of  and  too  small 
causes  the  same  type  of  computer  run  time  problems  as  dis- 
cussed above.  For  example,  a value  of  Tg=Tj=.001  results 
in  simulation  run  times  of  about  350  octal  seconds  for  a ten 
second  real  time  simulation.  Increasing  the  values  to 
Tg»TT*0.05  results  in  a simulation  run  time  of  45  octal  seconds. 

Since  the  recommended  range  of  values  of  Tg  and  for  the 
TAWDS  pilot  model  is  0.0  to  O.O5,  increasing  the  values  to  O.O5 
keeps  the  parameters  within  the  recommended  range. 

Two  other  considerations  for  the  use  of  the  EASY  programs 
are  in  order.  First,  the  setting  of  error  controls  for  the 
integrator  states  has  an  effect  on  the  accuracy  of  results 

J 

and  also  on  the  integration  efficiency.  The  default  value 

for  error  controls  is  0.1,  but  this  value  is  not  appropriate 

for  all  the  states.  The  values  of  the  state  outputs  vary 

widely!  for  instance,  the  altitude  state  output  might  be 

10,000  feet  while  the  output  of  a control  system  state  might 

be  0.1  degrees.  A more  appropriate  integrator  error  control  j 

i 

• i 

value  is  0.1  percent  of  the  estimated  maximum  value  of  the 
state.  The  second  consideration  is  the  type  of  integration 
routine  to  use.  Three  integration  methods  are  available  in 
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EASY I Fixed- Step  Euler,  Runga-Kutta,  and  Variable- Order  Gear. 
For  dynamic  simulation  of  nonlinear  systems,  the  Fixed  Step 
Euler  method  is  inappropriate  because  of  computation  errors. 
Runga-Kutta  and  Variable  Order  Gear  can  both  be  used.  Runga- 
Kutta  with  a maximum  step  size  of  0.2  seconds  was  used  in 
the  simulation  because  the  Runga-Kutta  method  proved  to  be 
more  efficient  due  to  the  nonlinearities  in  the  model. 


I 


j 

i 
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Appendix  B.  Parameter  Values  and  Aerodynamic  Data 


Appendix  B presents  a listing  of  the  constant  parameters 
and  tabular  aerodynamic  data  used  in  the  digital  simulation. 
This  data  was  provided  by  the  Air  Force  Flight  Dynamics 
Laboratory.  The  tabular  data  for  the  aerodynamic  coefficients 
was  obtained  by  curve  fitting  flight  test  data. 

The  weapon  delivery  constants  are  presented  in  the 
first  section.  The  second  section  presents  F-106b  airframe 
constants.  The  third  section  presents  the  F-106b  flight 
control  system  gain  schedules.  The  remaining  sections 
present  the  tabular  data  in  essentially  the  same  order  that 
the  data  are  used  in  the  digital  computer  program.  The 
EASY  program  acronyms  for  each  table  are  listed  in  paren- 
theses at  the  end  of  the  table  titles. 
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Weapon  Delivery  Constants 


1 


Bullet  time-of-f light  equation  constants* 

Ballistic  Coefficient,  = O.OO625 

Constants  assumed  to  represent  second  order  effects, 

c^  - 19.^ 

Cg  = 33^1 

Gun  angle  with  respect  to  aircraft  u axis, 

“ 2.5  degrees  (depressed  gun) 

Bullet  Muzzle  Velocity,  = 3350  ft/sec 

F-106b  Airframe  Constants 

Wing  Reference  Area,  S * 695  sq  ft 

Wingspan,  b = 38. 1 3 ft 

Mean  Aerodynamic  Chord,  c *=  23.755  ft 

Aircraft  Mass,  m = 924.7  slugs 

Moments  of  Inertia*  * 18,634  slug-ft^ 

ly  - 177,858  slug-ft^ 

- 191,236  slug-ft^ 

Product  of  Inertia*  I = 5539  slug-ft^ 

xz 

Distance  from  aircraft  eg  to  origin  of  the  thrust  vector* 
In  the  u direction,  = -11.8  ft 
In  the  w direction,  z^  » 0 
Tank  Buttock  Line,  (BL)^  = 9.37  ft 


Xcg  - -5.09  ft 

0 

1 

2cg  - -.76  ft 

^ref  ' 0 

'ref  * --^3  tX 

\ 

i 
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Flight  Control  System  Gain  Schedules 

Table  V.  Pitch  Rate  Feedback  Gain  (K^) 

Air  Density  Ratio  p/p^  0 ,05  .35  1 

Gain  -.4  -.4  -.2  -.1 

Pitch  Damper  Washout  Circuit  Time  Constant  ) 

q/ 

Air  Density  Ratio  p/p^  0 ,0?  .2  1 

Time  Constant  .5  ,5  .25  .05 


Table  VII » Dii'i'srential  Ele von- to- Rudder  Crossfeed  Gain  (K,  ) 

— — 5 . 


(Ib^ft^) 

0 

100 

150 

200 

300 

400 

500 

600 

800 

Gain  K, 

®A 

3 

3 

1.25 

.6 

-.2 

-.7 

-1.1 

-1.25 

-1. 

Table 

1 VIII 

. Roll 

Rate 

Feedback  Gain 

(Ib/ft^) 

0 

100 

150 

200 

300 

400 

500 

600 

800 

Gain  Kp 

.2  . 

155 

.13 

.11 

.08 

.065 

.05 

.042 

.033 
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Thrust  Required 


Table  IX.  Nominal  Thrust  Values  (lb) 


Altitude  (ft) 


Mach 

No. 

5000 

10000 

15000 

20000 

25000 

30000 

35000 

0.40 

9500 

8700 

7768 

6802 

5836 

5032 

4390 

0.88 

10300 

9350 

8515 

7565 

66I5 

5742 

4870 

0.90 

10400 

9450 

8584 

7624 

6664 

5794 

4924 

Lift  Coefficients 

Table  X.  Lift  Coefficient  at  a=0.  C,  (-CZO) 

~ ^ m . 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.0425 

-.0425 

-.0425 

-.0425 

-.0425 

0.50 

-.0485 

-.0586 

-.0487 

-.0488 

- . 0490 

0.80 

-.0560 

-.0573 

-.0580 

-.0584 

-.0588 

0.90 

-.0600 

-.0628 

-.0654 

-.0665 

-.0669 

1.00 

-.0575 

-.0637 

-.0686 

-.0719 

-.0738 

1.10 

-.0380 

-.0452 

-.0510 

-.0553 

-.0570 

1.20 

-.0262 

-.0324 

-.0375 

-.0414 

-.0426 

1.50 

-.0139 

-.0185 

-.0221 

-.0255 

-.0269 

1.80 

-.0127 

-.0166 

-.0203 

-.0230 

-.0243 

2.00 

-.0127 

-.0164 

-.0196 

-.0222 

-.0235 
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Table  XI.  Lift  Curve  Slope ^ C,  (-CZA) 

■^a 


Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

2.50c 

2.500 

2.500 

2.500 

2.500 

0.50 

2.560 

2.565 

2.570 

2.575 

2.580 

0.80 

2.720 

2.759 

2.780 

2.798 

2.802 

0.90 

2.860 

2.930 

2.970 

3.000 

3.018 

1.00 

3.079 

3.195 

3.290 

3.385 

3.420 

1.10 

2.910 

3.040 

3.170 

3.270 

3.310 

1.20 

2.680 

2.800 

2.940 

3.040 

3.095 

1.50 

2.145 

2.260 

2.370 

2.480 

2.525 

1.80 

1.800 

1.892 

1.989 

2.081 

2.125 

2.00 

1.625 

1.695 

1.800 

1.865 

1.912 

Table  XII.  Lift  Coefficient  Due  to  Elevator  Deflection. C.-  (-CZDE) 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.899 

.892 

.900 

.884 

.888 

0.50 

.830 

.882 

.903 

.921 

.929 

0.80 

.680 

.792 

.871 

.923 

.941 

0.90 

.536 

.679 

.799 

.870 

.897 

1.00 

.271 

.488 

.590 

.715 

.761 

1.10 

,200 

.309 

.431 

.540 

.598 

1.20 

.156 

.240 

.332 

.421 

.460 

1.50 

.080 

.120 

.179 

.231 

.262 

1.80 

.045 

.069 

.103 

.149 

.170 

2.00 

.034 

.050 

.082 

.111 

.130 
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Table 

XIII.  Lift 

Coefficient 

at  a*0  Due  to  Tanks,  i 

C,  (-CLOT) 

"°T 

Mach 

No. 

0.20 

0.50 

0.80  0.90 

1.0 

°T 

.000463 

.000289 

.000054  .000144 

.000226 

Mach 

No. 

1.1 

1.2 

1.5  1.8 

2.0 

°T 

.000235 

.000224 

.000210  .000194 

.000183 

Table  XIV.  Lift 

Curve  Slope 

Due  to  Tanks,  C,  (■ 

S 

-CLAT) 

Mach 

No. 

0.20 

0,50 

0.80  0.90 

1.0 

®T 

.000272 

.000253 

.000226  .000207 

.000189 

Mach 

No. 

1.1 

1.2 

1.5  1.8 

2.0 

\ 

“t 

.000180 

.000185 

.000180  .000180 

.000180 

Table  3^.  (CDRAG2) 

iJ4  ' 


Sg  (deg) 


Cl 

-15 

-10 

-5 

-2.5 

0 

-0.25 

.0395 

.0350 

.0312 

.0312 

.0340 

-0.20 

.0355 

.0300 

.0254 

.0254 

.0271 

-0.15 

.0320 

.0264 

.0210 

.0210 

.0219 

o 

ft 

. 

0 

1 

.0292 

.0237 

.0176 

.0176 

.0177 

-0.075 

.0281 

.0228 

.0165 

.0165 

.0165 

-0.05 

.0272 

.0219 

.0157 

.0154 

.0154 

-0.025 

.0266 

.0214 

.0152 

.0145 

.0145 

o 

• 

o 

.0262 

.0210 

.0149 

.0139 

.0139 

0.025 

.0262 

.0209 

.0147 

.0135 

.0135 

0.05 

.0265 

.0210 

.0149 

.0135 

.0134 

0.075 

.0280 

.0215 

.0152 

.0138 

.0136 

0.10 

.0304 

.0228 

.0159 

.0145 

.0140 

0.125 

.0333 

.0248 

.0168 

.0154 

.0148 

0.15 

.0375 

.0276 

.0183 

.0168 

.0157 

0.175 

.0425 

.0311 

.0206 

.0185 

.0170 

0.20 

.0480 

.0361 

.0236 

.0204 

.0189 

0.225 

.0535 

.0416 

.0280 

.0230 

.0208 

0.25 

.0600 

.0478 

.0330 

.0266 

.0235 

0.275 

.0660 

.0540 

.0384 

.0312 

.0266 

0.30 

.0730 

.0611 

.0440 

.0368 

.0305 

0.35 

.0874 

.0762 

.0565 

.0485 

.0410 

0.40 

.1028 

.0925 

.0706 

.0614 

.0535 

0.45 

.1185 

.1037 

.0854 

.0760 

.0684 

0.50 

.1340 

.1256 

.1003 

.0906 

.0840 
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Table 

XVI. 

^2 

(Mach>  0.9). 

(CDRAG3) 

(deg) 

^15 

-10 

-5 

-2.5 

0 

-0.25 

.0466 

.0400 

.0339 

.0329 

.0348 

o 

CM 

• 

0 

1 

.0421 

.0350 

.0285 

.0270 

.0280 

-0.15 

.0339 

.0305 

.0240 

.0221 

.0226 

-0.10 

.0354 

.0272 

.0205 

.0185 

.0185 

-0.075 

.0340 

.0259 

.0193 

.0171 

.0171 

-0.05 

.0331 

.0250 

.0181 

.0159 

.0159 

-0.025 

.0324 

.0242 

.0167 

.0154 

.0154 

o 

• 

o 

.0324 

.0238 

.0165 

.0149 

.0149 

0.025 

.0326 

.0238 

.0166 

.0147 

.0144 

0.05 

.0334 

.0242 

.0167 

.0149 

.0144 

0.075 

.0350 

.0248 

.0172 

.0151 

.0145 

0.10 

.0428 

.0303 

.0205 

.0175 

.0166 

0.125 

.0466 

.0336 

.0227 

.0189 

.0178 

0.15 

.0517 

.0383 

.0255 

.0210 

.0194 

0.175 

.0580 

.0430 

.0289 

.0238 

.0214 

0.20 

.0635 

.0482 

.0336 

.0275 

.0237 

0.225 

.0703 

.0536 

.0383 

.0316 

.0267 

0.25 

.0762 

.0600 

.0438 

.0365 

.0303 

0.275 

.0903 

.0731 

.0556 

.0476 

.0404 

0.30 

.1050 

.0880 

.0685 

.0607 

.0516 

0.35 

.1192 

.1035 

.0830 

.0755 

.0656 

0.40 

.1331 

.1188 

.0980 

.0910 

.0806 

0.45 

.1331 

.1188 

.0980 

.0910 

.0806 

0.50 

.1331 

.1188 

.0980 

.0910 

.0806 
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6g  (deg) 


-15 

-10 

-5 

-2.5 

0 

-0.25 

.0323 

.0294 

.0263 

.0275 

.0335 

-0.20 

.0285 

.0254 

.0223 

.0234 

.0264 

-0.15 

.0254 

.0216 

.0189 

.0196 

.0212 

-0.10 

.0230 

.0191 

.0160 

.0167 

.0175 

-0.075 

.0220 

.0182 

.0149 

.0156 

.0162 

-0.05 

.0215 

.0175 

.0139 

.0145 

.0150 

-0.025 

.0212 

, ''I72 

.0134 

.0139 

.0140 

0 

• 

0 

.0210 

.0171 

.0130 

.0134 

.0134 

0.025 

.0212 

.0171 

.0128 

.0132 

.0132 

0.05 

.0216 

.0174 

.0132 

.0132 

.0132 

0.075 

.0231 

.0181 

.0136 

.0136 

.0136 

0.10 

.0250 

.0191 

.0147 

.0141 

.0141 

0.125 

.0276 

.0213 

.0160 

.0151 

.0151 

0.15 

.0319 

.0241 

.0174 

.0162 

.0161 

0.175 

.0368 

.0294 

.0191 

.0176 

.0174 

0.20 

.0425 

.0334 

.0215 

.0194 

.0184 

0.225 

.0485 

.0393 

.0250 

.0214 

.0209 

0.25 

.0550 

.0454 

.0296 

.0245 

.0234 

0.275 

.0613 

.0520 

.0354 

.0290 

.0264 

0.30 

.0682 

.0590 

.0416 

.0345 

.0302 

0.35 

.0825 

.0735 

.0549 

.0466 

.0402 

0.40 

.0977 

.0890 

.0700 

.0615 

.0529 

0.45 

.1130 

.1045 

.0862 

.0785 

.0677 

0.50 

.1234 

.1196 

.1040 

.0975 

.0852 
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Pitchim 

t Moment 

Coefficients 

Table 

XVIII.  : 

Pitching  Moment 

Coefficient 

for  a=0. 

Cn,  (CMO) 

Altitude 

(ft) 

Mach 

0 

15000 

30000 

45000 

55000 

No. 

0.20 

.00410 

.00410 

.00410 

.00410 

.00410 

0.50 

.00487 

.00483 

.00487 

. 00492 

.00496 

0.80 

.00632 

.00636 

.00640 

.00644 

.00648 

0.90 

.00802 

.00820 

.00835 

.00840 

. 00845 

1.00 

.00972 

.01055 

.01099 

.01131 

.01162 

1.10 

.01088 

.01181 

.01251 

.01295 

.01316 

1.20 

.00805 

.00908 

.00955 

.00990 

.01009 

1.50 

.00245 

.00302 

.00359 

.00390 

.00161 

1.80 

.00028 

.00079 

.00121 

.00149 

.00161 

2.00 

-.00091 

-.00045 

.00005 

.00019 

.00042 

Table 

XIX.  Pitch 

Moment  Due 

to  Pitch 

Rate . C_ 

q 

(CMQ) 

Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.8100 

-.8800 

-.9000 

-.9300 

-.9600 

0.50 

-.5600 

-.6600 

-.7000 

-.7300 

-.7500 

0.80 

-.5400 

-.7100 

-.8000 

-.8600 

-.8800 

0.90 

-.5900 

-.8200 

-.9200 

-1.000 

-1.030 

1.00 

-.4900 

-.6800 

-.7500 

-.8400 

-.8900 

1.10 

-.3500 

-.4900 

-.5900 

-.6200 

-.6400 

1.20 

-.3600 

-.5000 

-.5900 

- . 6400 

-.6700 

1.50 

-.2900 

-.4100 

-.4800 

-.5300 

-.5400 

1.80 

-.2600 

-.3800 

-.4300 

-.4700 

-.4800 

2.00 

-.2600 

-.3800 

-.4200 

-.4500 

-.4700 
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Table 

XX.  Pitch 

Moment  Due  to 

Elevator 

Deflection 

C-  (CT.1DE) 

’"a 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.3350 

-.3300 

-.3370 

-.3320 

-.3322 

0.50 

-.3224 

-.3375 

-.3425 

-.3475 

- . 3482 

0.80 

-.3019 

-.3225 

-.3500 

-.3600 

-.3650 

0.90 

.02785 

-.3150 

-.3400 

-.3625 

-.3700 

1.00 

-.2287 

-.2837 

-.3285 

-.3607 

-.3770 

1.10 

-.1725 

-.2137 

-.2512 

-.2800 

-.2950 

1.20 

-.1362 

-.1650 

-.1950 

-.2212 

-.2312 

1.50 

-.0758 

-.0950 

-.1137 

-.1300 

-.1375 

1.80 

-.0475 

-.0600 

-.0750 

-.0875 

-.0938 

2.00 

-.0362 

-.0456 

-.0555 

-t0675 

-.0745 

Table 

XXI.  Pitch 

Moment  at 

a=0  Due  to 

Tanks  C_ 

(CMOT) 

°T. 

Mach 

No. 

0.20 

0.50 

0.80 

0.90 

1.0 

°T 

-.001142 

-.001151 

-.001169 

-.001151 

-.000941 

Mach 

No. 

1.0 

1.2 

1.5 

1.8 

2.0 

c.. 

°T 

-.000987 

-.001042 

-.001094 

-.001096 

-.001096 

Table 

XXII.  Pitch  Moment  Curve  Slope 

Due  to  Tanks 

(CMAT) 

Mach 

No. 

0.20 

0.50 

0.80 

0.90 

1.0 

“■"a 

®T 

.000158 

.000159 

.000159 

.000159 

.000160 

Mach 

No. 

1.0 

1.2 

1.5 

1.8 

2.0 

a 

®T 

.000160 

.000159 

.000161 

.000165 

.000166 

138 


Table  XXIII. 

Aircraft  Neutral  Point, 

N.P.  (TNP) 

Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.351 

-.351 

-.351 

-.351 

-.351 

0.50 

-.35^ 

-.354 

-.354 

-.354 

-.354 

0.80 

-.363 

-.365 

-.367 

-.369 

-.369 

0.90 

-.377 

-.382 

-.386 

-.388 

-.389 

1.00 

-.409 

-.420 

-.427 

-.434 

-.436 

1.10 

-.430 

-.442 

-.451 

-.460 

-.462 

1.20 

-.440 

-.453 

-.463 

-.470 

-.474 

1.50 

-.444 

-.456 

-.466 

-.475 

-.479 

1.80 

-.435 

-.446 

-.457 

-.465 

-.469 

2.00 

-.426 

-.439 

-.450 

-.458 

-.462 

Nonlinear  Side  Force  Coefficients 

Table  XXIV. 

Effect  of 

Tanks  on  C.. 

, C„  (CYBT) 
*^T 

Mach 

No. 

0.20 

0.50 

0.80 

0.90 

1.0 

-.00164 

-.00156 

-.00155 

-.00183 

-.00211 

Mach 

No. 

1.1 

1.2 

1.5 

1.8 

2.0 

-.00216 

-.00212 

-.00198 

-.00185 

-.00180 

Mach 

No. 

Table  XXV. 

0.20 

Effect  of  ' 

0.50 

Tanks  on  C . 

0.80 

C (CYAT) 

0.90 

1.0 

Cy 

S 

.001103 

.001021 

.000994 

.001111 

.001473 

Mach 

No. 

1.1 

1.2 

1.5 

1.8 

2.0 

Cv 

.001600 

.001609 

.001545 

.001491 

.001446 

1 
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(cm2) 


Table  XXVI. 


Second  Order  Effect  of  a on 


Altitude 

(ft) 

15000 

30000 

45000 

Mach 

No. 

0.20 

0.50 

0.80 

0.90 

1.00 

1.10 

1.20 

1.50 

1.80 

2.00 


Table 


Mach 

No. 

0.20 

0.50 

0.80 

0.90 

1,00 

1.10 

1.20 

1.50 

1.80 

2.00 


0 

-.000018 

-.000105 

-.000197 

-.000265 

-.000235 

-.000103 

-.000036 

-.000320 

-.000382 

-.000295 


0 

-.00483 

-.00881 

-.01400 

-.01559 

-.01409 

-.01138 

-.00939 

-.00505 

-.00197 

-.00075 


-.000018 

-.000105 

-.000197 

-.000265 

-.000235 

-.000103 

-.000036 

-.000330 

-.000438 

-.000378 


-.00483 

-.00881 

-.01420 

-.01569 

-.01420 

-.01171 

-.00963 

-.00509 

-.00191 

-.00063 


-.000018 

-.000105 

-.000197 

-.000265 

-.000235 

-.000103 

-.000036 

-.000334 

-.000461 

-.000418 


-.00483 

-.00881 

-.01426 

-.01579 

-.01441 

-.01206 

-.00975 

-.00512 

-.00185 

-.00050 


-.000018 

-.000105 

-.000197 

-.000265 

-.000235 

-.000103 

-.000036 

-.000341 

-.000467 

-.000431 


-.00483 

-.00881 

-.01432 

-.01589 

-.01460 

-.01220 

-.00987 

-.00516 

-.00177 

-.00038 


55000 

-.000018 

-.000105 

-.000197 

-.000265 

-.000235 

-.000103 

-.000036 

-.000337 

-.000463 

-.000429 


(CYBB) 


55000 

-.00483 

-.00881 

-.01438 

-.01599 

-.01478 

-.01234 

-.00999 

-.00519 

-.00171 

-.00025 


XXVII.  Second  Order  Effect  of  £ on  C • C 


Altitude  (ft) 

15000  30000 


45000 


PlPl 
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Table  XXVIII 

Effect  1 

of  a on  C . 

c„  (CYBAl 

Mach 

0.20 

0.50 

0.80 

0.90 

1.0 

No. 

.0051 

.0086 

.0142 

.0174 

.0198 

Mach 

1.1 

1.2 

1.5 

1.8 

2.0 

No. 

.0198 

.0182 

.0126 

.0084 

.0074 

Nonlinear  Rollins:  Moment  Coefficients 

Table  XXIX. 

Effect  of 

Tanks  on  C, 

, C,  (CLBT) 

X 

3 

Mach 

0.20 

0.50 

0.80 

0.90 

1.0 

No. 

.000101 

.000097 

.000096 

.000113 

.000131 

Mach 

1.1 

1.2 

1.5 

1.8 

2.0 

No. 

.000134 

.000131 

.000122 

.000114 

.000111 

Table 

XXX.  Second 

Order  Effect  of  a on 

(CLBA2) 

3 B 2 

Altitude  (ft) 

Mach 

0 

15000 

30000 

45000 

55000 

No. 

0.20 

.000156 

.000156 

.000156 

.000156 

.000156 

0.50 

.000152 

.000153 

.000153 

.000154 

.000154 

0.80 

.000144 

.000145 

.000146 

.000147 

.000148 

0.90 

.000113 

.000115 

.000117 

.000119 

.000121 

1.00 

.000089 

.000092 

.000094 

.000096 

.000096 

1.10 

.000138 

.000142 

.000146 

.000150 

.000154 

1.20 

.000180 

.000188 

.000196 

.000201 

.000204 

1.50 

.000226 

.000227 

.000227 

.000227 

.000228 

1.80 

.000235 

.000216 

.000201 

.000190 

.000186 

2.00 

.000204 

.000176 

.000153 

.000138 

.000129 

Table  XXXI . Effect  of  a on  C,  C,  ( CLBA ) 

Pf  "*3 


Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.00859 

-.00859 

-.00859 

-.00859 

-.00859 

0.50 

-.00870 

-.00872 

-.00874 

-.00877 

-.00879 

0.80 

-.00871 

-.00882 

-.00890 

-.00900 

-.00909 

0.90 

-.00858 

-.00877 

-.00892 

-.00901 

-.00911 

1.00 

-.00800 

-.00832 

-.00861 

-.00878 

-.00889 

1.10 

-.00695 

-.00726 

-.00755 

-.00777 

-.00786 

1.20 

-.00582 

-.00611 

-.00638 

-.00656 

-.00669 

1.50 

-,00322 

-.00340 

-.00357 

-.00369 

-.00379 

1.80 

-.00085 

-.00098 

-.00102 

-.00107 

-.00111 

2.00 

.00041 

.00044 

.00047 

.00050 

.00052 

Table  XXXII. 

Second  Order  Effect 

of  8 on  C- 

3. 

C (CLBB) 

3|el 

Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.00078 

-.00078 

-.00078 

-.00078 

-.00078 

0.50 

-.00143 

-.00143 

-.00144 

-.00144 

-.00144 

0.80 

-.00219 

-.00222 

-.00224 

-.00225 

-.00225 

0.90 

-.00231 

-.00236 

-.00239 

-.00242 

-.00243 

1.00 

-.00199 

-.00204 

-.00210 

-.00215 

-.00216 

1.10 

-.00152 

-.00161 

-.00165 

-.00170 

-.00172 

1.20 

-.00113 

-.00119 

-.00123 

-.00127 

-.00129 

1.50 

-.00021 

-.00022 

-.00023 

-.00024 

-.00024 

1.80 

.00041 

.00043 

.00044 

. 00045 

.00047 

2.00 

,00060 

.00064 

.00067 

.00069 

.00070 
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Table  XXXIII.  Second  Order  Effect  ^ a on  C.  « C.  (CLPA2) 

ip  Xp^2 

Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

450  DO 

55000 

0.20 

.000285 

.000286 

.000287 

.000288 

.000289 

0.50 

.000300 

.000307 

.000310 

.000312 

.000315 

0.80 

.000337 

.000359 

.000373 

.000382 

.000387 

0.90 

.000397 

.000430 

.000457 

.000459 

.000478 

1.00 

.000584 

.000644 

.000696 

.000740 

.000758 

1.10 

.000513 

.000576 

.000633 

.000677 

.000697 

1.20 

.000430 

.000486 

.000536 

.000577 

.000597 

1.50 

.000257 

.000295 

.000330 

.000360 

.000374 

1.80 

.000173 

.000199 

.000227 

.000248 

.000260 

2.00 

.000137 

.000161 

.000187 

.000204 

.000212 

Table  XXXIV 

. Effect 

of  a on  C.  • 
P 

C-  (CLPA) 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.00174 

-.00176 

-.00176 

-.00176 

^.00176 

0.50 

-.00176 

-.00179 

-.00183 

-.00184 

-.00185 

0.80 

-.00184 

-.00195 

-.00204 

-.00209 

-.00212 

0.90 

-.00200 

-.00216 

-.00229 

-.00238 

-.00242 

1.00 

-.00246 

-.00272 

-.00294 

-.00311 

-.OO3I8 

1.10 

-.00227 

-.00252 

-.00278 

-.00296 

-.00306 

1.20 

-.00202 

-.00228 

-.00253 

-.00273 

-.00281 

1.50 

-.00149 

-.00172 

-.00194 

-.00211 

-.00219 

1.80 

-.00117 

-.00136 

-.00155 

-.00171 

-.00178 

2.00 

-.00102 

-.00119 

-.00137 

-.00151 

-.00158 

L _J 


Table  XXXV. 

Effect 

of  a on  C,  , 
R 

C-  (CLRA) 

\ 

.Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.01455 

.01459 

.01464 

.01467 

.01470 

0.50 

.01480 

.01510 

.01540 

.01550 

.01560 

0.80 

.01615 

.01706 

.01775 

.01820 

.01835 

0.90 

.02058 

.02230 

.02310 

.02410 

.02485 

1.00 

.02135 

.02326 

.02520 

.02670 

.02730 

1.10 

.01429 

.01588 

.01780 

.01892 

.01944 

1.20 

.01100 

.01238 

.01370 

.01482 

.01558 

1.50 

.00755 

.00865 

.00970 

.01056 

.01097 

1.80 

.00609 

.00708 

.00802 

.00881 

.00920 

2.00 

.00548 

.00640 

.00730 

.00808 

.00843 

Nonlinear  Yawing  Moment  Coefficients 


Table  XXXVI.  Effect  of  Tanks  on  C_  • C (CNBT) 

Hq  tta  

P P J 


Mach 

No. 

0.20 

0.50 

0.80 

0.90 

l.U 

S 

Pt 

-.000300 

-.000282 

-.000264 

-.000219 

-.000180 

Mach 

No. 

1.1 

1.2 

1.5 

1.8 

2.0 

Prp 

-.000140 

-.000115 

-.000077 

-.000084 

-.000086 
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Table  XXXVII.  Second  Order  Effect  of  a on  C_  C_  (CNBAA) 

■* — na  np  ~ ' 


Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.000056 

-.000056 

-.000056 

-.000056 

-.000056 

0.50 

-.000066 

-.000066 

-.000066 

-.000066 

-.000066 

0.80 

-.000115 

-.000114 

-.000113 

-.000112 

-.000111 

0.90 

-.000147 

-.000146 

-.000145 

-.000144 

-.000143 

1.00 

-.000137 

-.000137 

-.000137 

-.000137 

-.000137 

1.10 

-.000058 

-.000058 

-.000058 

-.000058 

-.000058 

1.20 

-.000016 

-.000016 

-.000016 

-.000016 

-.000016 

1.50 

-.000242 

-.000242 

-.000242 

-.000242 

-.000234 

1.80 

-.000380 

-.000380 

-.000377 

-.000377 

-.000375 

2.00 

-.000435 

-.000399 

-.000385 

-.000380 

-.000385 

Table  XXXVIII.  Second  Order  Effect  of  B on  C_  C (CNBB) 

~ "b.  "s|b| 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.00235 

.00235 

.00235 

.00235 

.00235 

0.50 

.00420 

.00422 

.00425 

.00427 

.00430 

0.80 

.00650 

.00656 

.00659 

.00662 

.00665 

0.90 

.00665 

.00670 

.00680 

.00685 

.00690 

1.00 

.00580 

.00590 

.00605 

.00617 

.00621 

1.10 

.00420 

.00430 

. 00440 

. 00445 

.00450 

1.20 

.00260 

.00267 

.00274 

.00282 

.00290 

1.50 

-.00130 

-.00128 

-.00125 

-.00122 

-.00120 

1.80 

-.00370 

-.00375 

-.00380 

-.00385 

-.00390 

2.00 

-.00460 

-.00465 

-.00470 

-.00475 

-.00480 
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Table  XXXIX.  Effect  ^ a on  C_  . C (CNBA) 


Mach 

No. 

0,20 

0.50 

0.80 

0.90 

1.0 

Cn 

-.00303 

-.00259 

-.00195 

-.00165 

-.00310 

Mach 

1.1 

1.2 

1.5 

1.8 

2.0 

No. 

Cr, 

-.00293 

-.00311 

-.00318 

-.00267 

-.00239 

Table  XL. 

Effect  of 

a on  , 

*A 

C_  (CNDAA) 

"6, a 

A 

Mach 

No. 

0.20 

0.50 

0,80 

0.90 

1.0 

A 

.00074 

.00053 

-.00017 

-.00029 

.00094 

Mach 

No. 

1.1 

1.2 

1.5 

1.8 

2.0 

®n. 

.00179 

.00250 

.00241 

.00144 

.00108 

Table  XLI. 

Effect  of 

on 

E — n, 
*A 

f c 

"*A^E 

(CNDADE) 

Mach 

No. 

0,20 

0.50 

0.80 

0.90 

1.0 

^”6  6 
*A*E 

.00149 

.00165 

.00200 

.00236 

.00394 

Mach 

No. 

1.1 

1.2 

1.5 

1.8 

2.0 

^”6  6 
*A*E 

.00358 

.00294 

.00204 

.00147 

.00120 
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Linear  Lateral  Coefficients 


Table  XLII.  C (CYB) 
Altitude  (ft) 


Mach 

0 

15000 

30000 

45000 

55000 

No. 

0.20 

-.6060 

-.6060 

-.6060 

-.6060 

-.6060 

0.50 

-.6286 

-.6294 

-.6303 

-.6312 

-.6320 

0.80 

- . 6445 

- . 6494 

-.6510 

-.6529 

-.6548 

0.90 

-.6600 

-.6672 

-.6729 

-.6765 

-.6780 

1.00 

-.6461 

-.6590 

-.6700 

-.6780 

-.6819 

1.10 

-.6320 

-.6500 

-.6659 

-.6760 

-.6812 

1.20 

-.6188 

-.6390 

-.6562 

-.6675 

-.6729 

1.50 

-.5920 

-.6118 

-.6286 

- . 6420 

-.6471 

1.80 

-.5680 

-.5865 

-.6034 

-.6150 

-.6206 

2.00 

-.5550 

-.5719  -.5840 

Table  XLIII.  C 

Yp 

Altitude  (ft) 

-.6000 

(CYP) 

-.6058 

Mach 

0 

15000 

30000 

45000 

55000 

No. 

0.20 

-.00770 

-.00775 

-.00612 

-.00625 

-.00637 

0.50 

-.00475 

-.00550 

-.00600 

-.00625 

-.00637 

0.80 

-.00087 

-.00325 

-.00475 

-.00600 

-.00625 

0.90 

.00225 

-.00150 

-.00400 

-.00575 

-.00625 

1.00 

.00900 

.00300 

-.00150 

-.00475 

-.00600 

1.10 

.01 375 

.00575 

.00015 

-.00400 

-.00550 

1.20 

.01570 

.00745 

.00120 

-.00325 

-.00500 

1.50 

.01850 

.00950 

.00225 

• . 00250 

-.00475 

1.80 

.02010 

.01070 

.00325 

-.00175 

-.00400 

2.00 

..02P70 

.OHIO 

. .00350 

-.00137 

-.00375 

Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.6960 

.8960 

.8960 

.8960 

.8960 

0.50 

.9150 

.9160 

.9170 

.9190 

.9200 

0.80 

.9490 

.9560 

.9610 

.9640 

.9650 

0.90 

.9660 

.9650 

.9740 

.9800 

.9830 

1.00 

.7650 

.7920 

.8050 

.8140 

.8230 

1.10 

.7250 

.7480 

.7670 

.7800 

.7860 

1.20 

.7560 

.7800 

.8000 

.8170 

.8250 

1.50 

.6820 

.7050 

.7280 

.7460 

.7540 

1.80 

.5950 

.6140 

.6320 

.6470 

.6540 

2.00 

.5540 

.5710 

.5880 

.6020 

.6090 

■able  XLV.C  (CYDA) 
f>k 

Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.2040 

.2050 

.2060 

.2070 

.2080 

0.50 

.2220 

.2280 

.2320 

.2330 

.2330 

0.80 

.2200 

.2440 

.2590 

.2640 

.2680 

0.90 

.1960 

.2320 

.2630 

.2760 

.2840 

1.00 

.1360 

.2020 

.2530 

.2830 

.3030 

1.10 

.1020 

.1650 

.2240 

.2580 

.2810 

1.20 

.0600 

.1240 

.1740 

.2060 

.2280 

1.50 

.0140 

.0380 

.0580 

.0690 

.0760 

1.80 

.0070 

.0180 

.0300 

.0350 

.0390 

2.00 

.0070 

.0160 

.0260 

.0310 

.0340 
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Table  XL VI  r.  (CYDR) 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.1148 

.1148 

.1150 

.1150 

.1150 

0.50 

.1063 

.1100 

.1120 

.1131 

.1140 

0.80 

.0896 

.0987 

.1046 

.1080 

.1098 

0.90 

.0743 

.0852 

.0930 

.0973 

.0996 

1.00 

.0449 

.0572 

.0681 

.0745 

.0782 

1.10 

.0317 

.0419 

.0515 

.0588 

.0607 

1.20 

.0250 

.0340 

.0420 

.0480 

.0508 

1.50 

.0149 

.0218 

.0279 

.0330 

.0345 

1.80 

.0100 

.0159 

.0211 

.0256 

.0279 

2.00 

.0081 

.0140 

.0193 

.0238 

.0254 

Table  XLVII.n.  (CLP) 


JT 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.1650 

-.1670 

-.1666 

-.1670 

-.1670 

0.50 

-.1638 

-.1670 

-.1700 

-.1711 

-.1720 

0.80 

-.1601 

-.1693 

-.1762 

-.1809 

-.1826 

0.90 

-.1616 

-.1630 

-.1830 

-.1900 

-.1929 

1.00 

-.1670 

-.1834 

-.1975 

-.2088 

-.2139 

1.10 

-.1599 

-.1781 

-.1954 

-.2082 

-.2140 

1.50 

-.1403 

-.1612 

-.1802 

-.1957 

-.2030 

1.80 

-.1299 

-.1505 

-.1700 

-.1862 

-.1939 

2.00 

-.1234 

-.1438 

-.1639 

-.1801 

-.1875 
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Table  XLVIII 

■'r 

(CLR) 

Altitude 

(ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.138^ 

.1384 

.1384 

.1384 

.1384 

0.50 

.1333 

.1333 

.1333 

.1333 

.1333 

0.80 

.1172 

.1177 

.1178 

.1179 

.1180 

0.90 

.1119 

.1123 

.1123 

.1124 

.1125 

1.00 

.0885 

.0893 

.0901 

.0909 

.0917 

1.10 

.1112 

.1159 

.1198 

.1224 

.1234 

1.20 

.1084 

.1126 

.1163 

.1194 

.1203 

1.50 

.0933 

.0965 

.0995 

.1014 

.1022 

1.80 

.0832 

.0862 

.0884 

.0901 

.0908 

2.00 

.0795 

.0822 

.0845 

.0851 

.0877 

Table  XLIX.  C. 

P 

(CLB) 

Altitude  (ft) 

Mach 

No. 

0 

15000  30000 

45000 

55000 

0.20 

-.0280 

-.0280  -.0280 

-.0280 

-.0280 

0.50 

-.0279 

-.0278  -.0276 

-.0275 

-.0274 

0.80 

-.0202 

-.0196  -.0187 

-.0183 

-.0180 

0.90 

-.0188 

-.0178  -.0156 

-.0165 

-.0160 

1.00 

-.0291 

-.0284  -.0279 

-.0272 

-.0266 

1.10 

-.0524 

-.0544  -.0559 

-.0571 

-.0580 

1.20 

-,0588 

-.0618  -.0645 

-.0656 

-.0663 

1.50 

-.0562 

-.0583  -.0610 

-.0624 

-.0632 

1.80 

-.0483 

-.0507  -.0528 

-.0539 

-.0542 

2.00 

-.0404 

-.0425  -.0444 

-.0450 

-.0454 
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Table  L.p  (CLDA) 

* c 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.1140 

-.1180 

-.1200 

-.1200 

-.1200 

0.50 

-.1080 

-.1160 

-.1200 

-.1230 

-.1250 

0.80 

-.0830 

-.0980 

-.1120 

-.1200 

-.1250 

0.90 

-.0590 

-.0810 

-.0980 

-.1080 

-.1140 

1.00 

-.0190 

-.0400 

-.0600 

-.0740 

-.0830 

1.10 

-.0090 

-.0230 

-.0390 

-.0530 

-.0620 

1.20 

-.0070 

-.0190 

-.0320 

-.0440 

-.0510 

1.50 

-.0040 

-.0120 

-.0230 

-.0330 

-.0390 

1.80 

-.0010 

-.0080 

-.0160 

- . 0240 

-.0290 

2.00 

-.0010 

-.0050 

-.0120 

-.0180 

-.0220 

Table  LI. 

C,  (CLDR) 

K 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.0209 

.0209 

.0209 

.0206 

.0206 

0.50 

.0190 

.0197 

.0202 

.0206 

.0208 

0.80 

.0156 

.0175 

.0188 

.0196 

.0200 

0.90 

.0125 

.0148 

.0165 

.0175 

.0179 

1.00 

.0066 

.0094 

.0116 

.0133 

.0140 

1.10 

.0046 

.0068 

.0089 

.0102 

.0109 

1.20 

.0034 

.0054 

.0072 

.0086 

.0091 

1.50 

.0018 

.0032 

.0046 

.0056 

.0062 

1.80 

.0009 

.0022 

.0034 

.0046 

.0049 

2.00 

.0006 

.0019 

.0030 

.0040 

.0045 

L 
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Table  LII.  C (CNB) 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.1090 

.1090 

.1090 

.1090 

.1090 

0.50 

.1148 

.1151 

.1154 

.1157 

.1160 

0.80 

.1210 

.1235 

.1250 

.1264 

.1270 

0.90 

.1237 

.1268 

.1299 

.1324 

.1325 

1.00 

.1408 

.1477 

.1538 

.1570 

.1590 

1.10 

.1490 

.1580 

.1660 

.1720 

.1740 

1.20 

.1428 

.1526 

.1617 

.1685 

.1704 

1.50 

.1164 

.1262 

.1350 

.1405 

.1437 

1.80 

.0949 

.1036 

.1112 

.1173 

.1200 

2.00 

.0834 

.0915 

.0994 

.1049 

.1074 

Table  LIII 

(CNDA) 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.0720 

-.0720 

-.0730 

-.0730 

-.0730 

0.50 

-.0730 

-.0750 

-.0760 

-.0770 

-.0780 

0.80 

-.0760 

-.0830 

-.0880 

-.0900 

-.0920 

0.90 

-.0710 

-.0870 

-.0970 

-.1020 

-.1050 

1.00 

-.0560 

-.0850 

-.1100 

-.1240 

-.1340 

1.10 

-.0400 

-.0730 

-.0990 

-.1150 

-.1260 

1.20 

-.0280 

-.0540 

-.0780 

-.0920 

-.1020 

1.50 

-.0070 

-.0170 

-.0250 

-.0320 

-.0360 

1.80 

-.0030 

-.0090 

-.0140 

-.0170 

-.0190 

2.00 

-.0040 

-.0080 

-.0120 

-.0140 

-.0160 

Table  LIV.  (CNDR) 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.0570 

-.0570 

-.0570 

-.0570 

-.0570 

0.50 

-.0535 

-.0550 

-.0560 

-.0565 

-.0570 

0.80 

-.0463 

-.0505 

-.053^ 

-.0552 

-.0556 

0.90 

-.0405 

-.0457 

-.0495 

-.0517 

-.0527 

1.00 

-.0278 

-.0347 

-.0400 

-.0436 

-.0453 

1.10 

-.0203 

-.0259 

-.0307 

-.0346 

-.0359 

1.20 

-.0160 

-.0210 

-.0250 

-.0288 

-.0298 

1.50 

-.0100 

-.0137 

-.0171 

-.0196 

-.0209 

1.80 

-.0076 

-.0104 

-.0134 

-.0158 

-.0170 

2.00 

-.0063 

-.0091 

-.0120 

-.0145 

-.0155 

Table  LV. 

C (CNP) 

Hp 

Altitude  (ft) 

Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

.0073 

.0073 

.0073 

.0074 

.0074 

0.50 

.0068 

.0070 

.0073 

.0075 

.0075 

0.80 

.0055 

.0066 

.0074 

.0079 

.0081 

0.90 

.0045 

.0063 

.0076 

.0082 

.0086 

1.00 

.0021 

.0052 

.0075 

.0091 

.0096 

1.10 

.0004 

.0041 

.0071 

.0090 

.0098 

1.20 

-.0006 

.0034 

.0066 

.0088 

.0097 

1.50 

-.0019 

.0024 

.0058 

.0082 

.0092 

1.80 

-.0028 

.0016 

.0051 

.0075 

.0085 

2.00 

-.0032 

.0012 

.0046 

.0071 

.0082 

Table  LVI.  (CNR) 

R 


Altitude  (ft) 


Mach 

No. 

0 

15000 

30000 

45000 

55000 

0.20 

-.378 

-.378 

-.378 

-.378 

-.378 

0.50 

-.400 

-.402 

-.402 

-.403 

-.403 

0.80 

-.430 

-.^33 

-.436 

-.438 

-.^39 

0.90 

-.443 

-.448 

-.453 

-.456 

-.458 

1.00 

-.374 

-.384 

-.392 

-.395 

-.398 

1.10 

“.356 

-.367 

-.377 

-.383 

-.387 

1.20 

-.369 

-.382 

-.394 

-.402 

-.406 

1.50 

-.338 

-.3^9 

-.361 

-.368 

-.372 

1.80 

-.299 

-.309 

-.318 

-.324 

-.328 

2.00 

-.281 

-.290 

-.297 

-.304 

-.307 
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